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Abstract 

(D : 

'T^ , In the high-energy Umit, we compute the gauge-invariant three-parton forward 

.5^ I clusters, which in the BFKL theory constitute the tree parts of the NNLO impact 

^ ' factors. In the triple collinear limit, we obtain the polarized double-splitting func- 

d ' tions. For the unpolarized and the spin-correlated double-splitting functions, our 

results agree with the ones obtained by Campbell-Glover and Catani-Grazzini, re- 
spectively. In addition, we compute the four-gluon forward cluster, which in the 
BFKL theory forms the tree part of the NNNLO gluonic impact factor. In the 
quadruple collinear limit we obtain the unpolarized triple-splitting functions, while 
in the limit of a three-parton central cluster we derive the Lipatov vertex for the 
production of three gluons, relevant for the calculation of a BFKL ladder at NNLL 
accuracy. Finally, motivated by the reorganization of the color in the high-energy 
limit, we introduce a color decomposition of the purely gluonic tree amplitudes in 
terms of the linearly independent subamplitudes only. 



1 Introduction 



QCD calculations of multijet rates beyond the leading order (LO) in the strong coupling 
constant are generally quite involved. However, in recent years it has become clear 
how to construct general-purpose algorithms for the calculation of multijet rates at next- 
to-leading order (NLO) accuracy [|], 0, H, ||]. The crucial point is to organise the 
cancellation of the infrared (i.e. coUinear and soft) singularities in a universal, i.e. process- 
independent, way. The universal pieces in a NLO calculation are given by the tree-level 
splitting and eikonal 0, |[ functions, and by the universal structure of the poles of the 
one-loop amplitudes ^, |^. 

Eventually, the same procedure will permit the construction of general-purpose al- 
gorithms at next-to-next-to-leading order (NNLO) accuracy. It is mandatory then to 
fully investigate the infrared structure of the phase space at NNLO. The universal pieces 
needed to organise the cancellation of the infrared singularities are given by the tree-level 



double-splitting pOlO, double-eikonal [H, 12] and splitting-eikonal |T0, O functions, by 



the one-loop splitting [|I3|, and eikonal [|13[ functions, and by the universal structure 
of the poles of the two-loop amplitudes |T^. 

Another outstanding issue in QCD, at first sight unrelated to the topics discussed 
above, is the calculation of the higher-order corrections to the BFKL equation [|T^, |T^. 
In scattering processes characterised by two large and disparate scales, like s, the squared 
parton center-of-mass energy, and t, a typical momentum transfer, the BFKL equation 
resums the large logarithms of type ln(s/t). The LO term of the resummation requires 
gluon exchange in the cross channel, which for a given scattering occurs at 0{a1). The 
corresponding QCD amplitude factorizes then into a gauge-invariant effective amplitude 
formed by two scattering centers, the LO impact factors, connected by the gluon ex- 
changed in the cross channel. The LO impact factors are characteristic of the scattering 
process at hand. The BFKL equation resums then the universal leading-logarithmic (LL) 
corrections, of ln"(s/t)), to the gluon exchange in the cross channel. The building 



blocks of the BFKL resummation are the Lipatov vertex ||18[, i.e. the effective gauge- 



invariant emission of a gluon along the gluon ladder in the cross channel, and the gluon 
reggeization i.e. the LL part of the one- loop corrections to the gluon exchange in the 
cross channel. 

The accuracy of the BFKL equation is improved by computing the next-to-leading 
logarithmic (NLL) corrections |jl9|, |20|, i.e. the corrections of ln"~^(s/t)), to the 



gluon exchange in the cross channel. In order to do that, the universal building blocks 
of the BFKL ladder must be computed to NLL accuracy. These are given by the tree 
corrections to the Lipatov vertex, i.e. the emission of two gluons |2^, ^ or of a qq 
pair [^, 24 1 along the gluon ladder, by the one- loop corrections to the Lipatov vertex |25, 
6|, and finally by the NLL gluon reggeization [^, i.e. the NLL part of the two- 



loop corrections to the gluon ladder. However, to compute jet production rates at NLL 
accuracy, the impact factors must be computed at NLO [29, |3^. For jet production at 
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large rapidity intervals, they are given by the one-loop corrections to the LO impact 
factors, and by the tree corrections pTl |2^, E!], E^, i.e. the emission of two partons in 



the forward-rapidity region. In the collinear or soft limits, the latter reduce to the tree 
splitting or eikonal functions [0^] . 

To further improve the accuracy of the BFKL ladder one needs to compute the next-to- 
next-to-leading logarithmic (NNLL) corrections, i.e. the corrections of ln"~^(s/t)), 
to the gluon ladder. At present it is not known whether such corrections can be resummed. 
If that is the case, the universal building blocks of a BFKL ladder at NNLL would be: the 
emission of three partons along the gluon ladder, the one-loop corrections to the emission 
of two partons along the ladder, the two-loop corrections to the Lipatov vertex, and the 
gluon reggeization at NNLL accuracy. None of them is known at present. In this paper 
we compute the gluonic NNLO Lipatov vertex, i.e. the emission of three gluons along the 
ladder. 

In addition, to compute jet production rates at NNLL accuracy, the BFKL ladder 
should be supplemented by impact factors at NNLO. They are not known either. In 
this paper we compute their tree components, i.e. the emission of three partons in the 
forward-rapidity region. By taking then the triple collinear limit of the tree NNLO impact 
factors, we obtain the polarized double-splitting functions. Summing over the parton po- 
larizations, we obtain the unpolarized and the spin-correlated double-splitting functions, 
previously computed in Ref. [|T^] and [|TT|, respectively, in the conventional dimensional 
regularization (CDR) scheme. Since we sum over two helicity states of the external par- 



tons, as it is done in the dimensional reduction (DR) scheme 0], our results agree 



with the ones in the CDR scheme by setting there the dimensional regularization scheme 
(RS) parameter e = 0. 

For a scattering with production of m partons, we define the n-parton cluster, with 
m > n, as the set of n final-state partons where the distance in rapidity between any two 
partons in the cluster is much smaller than the rapidity distance between a parton inside 
the cluster and a parton outside. In the BFKL theory, {n + l)-parton forward clusters 
provide the tree parts of N"LO impact factors, while [n + l)-parton central clusters 
provide the tree parts of the N^LO Lipatov vertex, n-parton clusters were given also a 
field-theoretical basis in terms of an effective action describing the interaction between 
physical partons grouped into gauge-invariant clusters and the gluons exchanged in the 
cross channel In addition to computing the three-parton forward clusters and the 
three-gluon central cluster mentioned above, we compute the four-gluon forward cluster, 
i.e. the purely gluonic tree part of the NNNLO impact factor. By taking then the 
quadruple collinear limit, we obtain the polarized triple-splitting functions. They could 
be used in a gauge-invariant evaluation of the Altarelli-Parisi evolution at three loops |^ 



The outline of the paper is: in Section ^ we review the standard color decompositions 
of the n-parton tree amplitudes, and we present a color decomposition of the gluon am- 
plitudes in terms of the linearly independent subamplitudes only. In Section ^ we review 
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the elastic scattering of two partons in the high-energy hmit, which allows for the extrac- 
tion of the LO impact factors. In Section ^ we review the amplitudes for the production 
of three partons, with a gauge-invariant two-parton forward cluster; from these, we can 
extract the tree parts of the NLO impact factors; by taking the coUinear limit, we obtain 
the LO splitting functions. In Section |^ we compute the amplitudes for the production 
of four partons, with a three-parton forward cluster; then we extract the tree parts of the 
NNLO impact factors, and by taking the triple coUinear limit we obtain the polarized and 
unpolarized double-splitting functions. In Section we compute the amplitude for the 
production of five gluons, with a four-gluon forward cluster. We extract the tree part of 
the gluonic NNNLO impact factor, and by taking the quadruple coUinear limit we obtain 
the polarized triple-splitting functions. In addition, by taking the limit in which three glu- 
ons are emitted in the central-rapidity region, we obtain the gauge-invariant three-gluon 
central cluster, i.e. the tree part of the NNLO Lipatov vertex. In Section |^ we draw our 
conclusions. 



2 Tree Amplitudes 

In this section we review the color decomposition of purely gluonic and quark-gluon tree 
amplitudes. For the purely gluonic tree amphtudes, we introduce a color decomposition 
in terms of the linearly independent subamplitudes, Eq. (|2.9|). 

2.1 Gluon amplitudes 



For an amplitude with n gluons the usual color decomposition at tree level reads [^, |3^, 

1113,11, 

iA{gi,...,gn) = tg""-' Yl tr(A'^^ ■■ ■A'^")A((?.„ ... (2.1) 

where Sn/'^n are the non-cyclic permutations of n elements. The dependence on the 
particle helicities and momenta in the subamplitude, and on the gluon colors in the trace, 
is implicit in labelling each leg with the index i. Helicities and momenta are defined as if 
all particles were outgoing. 



The gauge invariant subamplitudes A satisfy the relations |40, 41 1, proven for arbitrary 



n in Ref . E2 



y4(l, 2, . . . , n — 1, n) = A(n, 1, 2, . . . , n — 1) cychcity 
A(l,2,...,n) = (-l)"A(n,...,2,l) reflection (2.2) 

A(l, 2, 3, . . . , n) + A(2, 1, . . . , n) + . . . + A(2, 3, . . . , 1, n) = dual Ward identity 

The above relations are sufficient to show that, for n < Q the number of independent 
subamplitudes can be reduced from (n — 1)! to {n — 2)\. For n > 7 it is still possible to 
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introduce a basis of (ra — 2)! elements by using Kleiss-Kuijf 's relation [^| 

x^,..., Xp, 2, 1/1, ... , y^) = {-If ^ A{1, 2, {a}{(3}) (2.3) 

crGOP{a}{/3} 

where e {a} = {a;^, Xp^i, . . . ,xi} , Pi G = {yi, . . . , j/g} and 0P{q;}{/5} is the set 
of permutations of the {n — 2) objects {xi, . . . ,Xp,yi, . . . , yq} that preserve the ordering 
of the Oi within {a} and of the Pi within {P}, while allowing for all possible relative 
orderings of the with respect to the Pi. The above relation has been checked up to 
n = 8 in Ref . , and proven for arbitrary n in Ref . |Q . Accordingly, the expression for 
the summed amplitude squared can be written as 



n] 



{n-iy. 

CijAiA* 

Cum J2 

(n-2)! 

cijA-iA* . 



|A(l,a2,...,a„)|2 + 



1 



where Cj, in Eq 



IS 



= (Z)'^-' E tK^.(A'S . . . , A'^"))[tr(P,(A'^\ . . . , A'^"))]* 



colors 



with Pj the i permutation in Snf^n- In Eq- (P-5|), the coefficient Cn{Nc) is 



C„(iV, 



{Nl - 1) 



(2.4) 
(2.5) 

(2.6) 



(2.7) 



(2.8) 



The first term in Eq. ( |2.5| ) constitutes the Leading Color Approximation (LCA). Up to 
n = 5, the l/N"^ corrections in Eq. ( p.5|) vanish and LCA is exact. The reduced color 
matrix Cij in Eq. (|2.6|) , has been obtained from Cjj applying the linear transformations 
of Eq. (|2.3|) , thus the labels z, j in Eq. (|2.6|) run only on the permutations of the linearly 
independent subamplitudes. 



Motivated by the reorganization of the color in the high-energy limit ^ 
and using Eqs. ( p.2|) and ( |2.3| ) we rewrite Eq. ( p.l|) as 



^(5-1, . . . ,5^^ 



n-2 



J J ^\yii 9(721 ■■■ 1 9(7,1-11 Qn 



-.n-2 



J2 (P-^ . . .p'^"-i),^,„A((?i, (7^,, . . .,g^^_„g„ 



(2.9) 



(7&Sn 



where {F"-)bc = z/''"'^. We have checked Eq. ( |2^) up to = 7. Eq. ( |2.9|) enjoys several 
remarkable properties. Firstly, it shows explicitly which is the color decomposition that 
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allows us to write the full amplitude iA in terms of the {n — 2)\ linearly independent 
subamplitudes only. In the following we shall refer to it as to a color ladder. Hence the 
color matrix obtained squaring Eq. ( |2.9| ) yields directly the Cij matrix in Eq. ([^.61). We 
have checked it against the explicit results of Ref. |4^, up to n = 5. Moreover, it is 



quite suggestive to note the formal correspondence with the amplitudes with a quark- 
antiquark pair and {n — 2) gluons, Eq. where the only difference between the two[| 

is the appropriate representation for the color matrices, namely the adjoint for the ra-gluon 
amplitude and the fundamental for the one with the qq pair. Finally, the most relevant 
applications of Eq. (|2.9|) for this work are to the study of the multi-gluon amplitudes in the 



high-energy limit. As discussed in the following, the color ladder naturally arises [^, 57 
in the configurations where the gluons are strongly ordered in rapidity, i.e. in the multi- 
Regge kinematics. Indeed in the strong-rapidity ordering only the subamplitude with the 
corresponding order in the color coefficient contributes to Eq. (|2.9|) . At NLO, where the 
strong ordering is relaxed for two adjacent gluons, the leading subamplitudes are the two 
which differ just by the exchange of the gluon labels in the color ladders 0]. As we shall 
see this result generalizes at NNLO and beyond. Nonetheless, in the following we have 
chosen to derive our results starting from Eq. ( |2.1| ) instead of using directly Eq. ( p.9| ). The 
former, though more laborious, shows explicitly how the color traces must be recombined 
to obtain the color ladder and, more importantly, allows us to find the relations necessary 
to prove the factorization in the multi-coUinear limits. 

For the maximally helicity-violating configurations, (— , — , -f, +), in Eq. (|2.1| ) or 
Eq. ( |2.9|) , there is only one independent color/helicity subamplitude, the Parke- Taylor 
(PT) subamplitude 

A{g,, g^) = 2-/2— , (2.10) 

(1 2) ■ ■ ■ {[n — 1) n) {n 1) 

where the i*^ and the j^'^ gluons have negative helicity. All other color/helicity amplitudes 
can be obtained by relabelling and by use of refiection symmetry, Eq. ( p2.2| ), and parity 
inversion. Parity inversion fiips the helicities of all particles, and it is accomplished by the 
substitution (ij) ^ [j i]. Subamplitudes of non-PT type, i.e. with three or more gluons 
of — helicity have a more complicated structure. 



2.2 Quark-gluon amplitudes 

For an amplitude with two quarks and {n — 2) gluons the color decomposition at tree-level 

is 13,11, n, 

^A{q,q■,g,,...,9in-2))=^9''^' (A'^^ .• A'^"-)/ g; , (2.11) 

a£Sn-2 

*The factor 1/2 in front of Eq. (|2.9| ) is due to our choice for the normahzation of the fundamental 
representation matrices, i.e. tr(A°A'') — S°'''/2. 
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where Sn-2 is the permutation group on n — 2 elements. 

For the maximally helicity-violating configurations, (—, —,+,..., +), there is one in- 
dependent color/hehcity subamphtude, the Parke- Taylor (PT) subamplitude 



\ 9(n-2)/2 
,g(n-2)) =2^ " - 



(2.12) 



(gg)(gl)---((n-2) g) ' 

where gluon gi has negative helicity. Helicity is conserved along the massless-fermion line. 
All other color /helicity amplitudes can be obtained by relabelling and by use of parity 
inversion, refiection symmetry and charge conjugation. In performing parity inversion, 
there is a factor of —1 for each pair of quarks participating in the amplitude. Refiec- 
tion symmetry is like in Eq. ( |2.2| ), for gluons and/or quarks alike. Charge conjugation 
swaps quarks and antiquarks without inverting helicities. In particular, using refiection 
symmetry and charge conjugation on Eq. ( p.l2| ) we obtain 

{qif{qi) 



A(g-,g+;^i,...,^(„_2)) = 2("-^)/2- 



(2.13) 



(gg)(gl)---((n-2) q) 

where gluon Qi has negative helicity. 

For an amplitude with four quarks and (n — 4) gluons the color decomposition at 
tree-level is 



n-4 



iA{qi,qi]q2,q2]gi, 



,9{n-4:)) 



^9' 



(2.14) 



X 



(A-\..A'^'=)^.;^ iXP^...y')/^^ Aiq,,qr,q2,q2;9a„...,9.,;9pr 



•9p, 



1 



(A'^\..A'^^),:^ (A^^..A^O,? B{q-,,qi;q-2,q2;9.^,-,9a,;9p^ 



'9 PI 



with /c -f / = n — 4, and where we suppose that the two quark pairs have distinct fiavor. 
The sums are over the partitions of [n — 4) gluons between the two quark lines, and over 
the permutations of the gluons within each partition. For = or / = 0, the color strings 
reduce to Kronecker delta's. For identical quarks, we must subtract from Eq. ( p.l4| ) the 
same term with the exchange of the quarks (gi ^ ^2)- 

For the maximally helicity-violating configurations, (— , — , +, -|-), with like-helicity for 
all of the gluons, the A and B subamplitudes factorize into distinct contributions for the 
two quark antennae |3^, Bol El|]. However, as we shall see in Sect. p.4|, we need the 



helicity configurations with two gluons of opposite helicity. For these the above mentioned 
factorization does not occur. 



3 The Leading Impact Factors 

We consider the elastic scattering of two partons of momenta pa and pb into two partons of 
momenta pa' and py, in the high-energy limit, s ^ Firstly, we consider the amplitude 
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Pa a l^a 



Pa' a Va 



2 



(a) 




Figure 1: (a) Amplitude ioi g g ^ g g scattering and (6), (c) for qg ^ qg scattering. We 
label the external lines with momentum, color and helicity, and the internal lines with 
momentum and color. 



for gluon-gluon scattering (Fig. |I]a). Using Eqs. ( p.l| ), (pT2|), or Eq. (p.9|), and Eq. ( p.lO| ), 
and Appendix we obtain 

A''^''{p:\py I Pir^pi") = 2^ [wr''c^''{pT;py)] \ [^9f''''c''''{pl''■,p7)] > (s.i) 



with q = Pb' + Pb and t 
gluon, are 



|g_Lp. The LO impact factors g* g g, with (?* an off-shell 



C'^iPa^Pt'] 



C'-'^iPb^Pt') 



Pb'. 



Pb'. 



(3.2) 



They conserve helicity along the on-shell gluon line and transform under parity into their 
complex conjugates, 

[C^';9({r})]* = C3'»({A;~"}) . (3.3) 

In Eq. ( p . 1| ) four helicity configurations are leading, two for each impact factor f|. The 
helicity-flip impact factor C^''^{p'^]p''^) is sub leading in the high-energy limit. 

From Eqs. ( [^.11D -( P.12D , we obtain the quark-gluon qg ^ qg scattering amplitude in 
the high-energy limit [ PT[] , 

A'''-''\p:%py I = 2s [9K'-aC^'''iPa''^';py)] j [^9f'''c^'\Pb';p:n] m 
A''^''{p:%py I py^Pi") = 2s [tgr''c^''{p:'';py)] ] [9K'-bC'-''{p7''';p7)] (3.5) 

where we have labelled the incoming quarks as outgoing antiquarks with negative mo- 
mentum, e.g. the antiquark is pa in Eq. p.4| ) (Fig. and pb in Eq. ( |3.5| ) (Fig. |I|c). The 
LO impact factors g* q ^ q are. 



c'-''iPa;py 



c'"^{Pb;pt') = ^ 



Pb'. 



Pb'. 



1/2 



(3.6) 



'''All throughout this paper, we shall always write only half of the helicity configurations contributing 
to an impact factor, the other half being obtained by parity. 
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Under parity, the functions ( |3.(j| ) transform as 

= SC^'^iik-'}) with S = -sign(g%°) 



(3.7) 



and in general an impact factor acquires a coefficient 5* for each pair of quarks (see Sect. ^. 
Analogously, the antiquark-gluon qg ^ qg amplitude is 



A 



''^'\p:\py I p7,pI') = 2s [gK-a,C'^'%p-''^';py)] ^ igf'''C^-\pl\p7) (3.8) 



A- 



gq^gq 



\Pa ' Pa' I Pb' ' Pb 



2s 



t 



igr''C^''{p7;p7) - gKi'C'-''{pr';pir) (3.9) 



where the antiquark is Pa' in Eq. ( p.8|) and in Eq. (|3.9| ), and the LO impact factors 
g* q ^ q are. 



C''-'%p-;p^,)=z; C™>,-,p+) = -^ 



Pb'± 
Pb'± 



1/2 



(3.10) 



In the amplitudes (|3.1| ), ( p.4|) , ( p.5|) , (p.8|) , (|3l9|) , the leading contributions from all the 
Feynman diagrams have been included. However, the amplitudes have the effective form 
of a gluon exchange in the t channel (Fig. p, and differ only for the relative color strength 



in the production vertices |^8[. This allows us to replace an incoming gluon with a quark, 
for instance on the upper line, via the simple substitution 



^g^^c^■''ip:-,py) --gK'-aC'-''{p7'^':p7 



(3.11) 



and similar ones for an antiquark and/or for the lower line. For example, the quark-quark 
qq qq scattering amplitude in the high-energy limit is 



A'^-''{p:'^,py I p7,pI'') = 2s [gK,,C^''iP~''^';py)] - 



9KiC'-'%pP';p7: 



. (3.12) 



4 The Next-to-leading Impact Factors 

Let three partons be produced with momenta ki, k2 and pb> in the scattering between two 
partons of momenta pa and pb, and to be specific, we shall take partons ki and /c2 in the 
forward-rapidity region of parton pa, the analysis for ki and k2 in the forward-rapidity 
region of pf, being similar. Parametrizing the momenta as in Eq. ( |A.1| ), we have, 

Z/i ^ 1/2 > Z/fe' ; \ki±\^\k2±\^\pb'±\. (4.1) 



4.1 The NLO impact factor g g* g g 

We consider the amplitude for the scattering gg ^ ggg (Fig. |^a). Only PT subam- 
plitudes contribute, thus using Eqs. (p?T|), (|2.2| ) and ( |2.10|) , and Appendix y, we ob- 
tain 0, m 

A^'-^'{p':%ki\k''^-\pyy,^) 
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Pa a Va 



ki di V{ 



o 
o 



^Tnr^5i5ir5wciro7rQip ^7^^w^^l5w^^o75l5l5^p ^rnfwsiswfruTrsip 



Pbb i^b , . Pb' b' 

[a) 



Figure 2: Amplitudes for the production of three partons, with partons ki and k2 in the 
forward-rapidity region of parton pa- 



(T2 ^ 
0"2 , 



(4.2) 



xtr (X'X'^^^X^^^X^'X^ - X'^X^-^X'^^^X^X^' + A'^A^'A^A'^'^^A'^-i - A'^A^A^'A'^'^^ A"^"^!) 
+B3-^39^pUa. kl'[\kZ^)iT (X'^X'^^^X^'X^X'^^^ - X'^X'^^^X^X^'X'^^^^ 



with the sum over the permutations of the two gluons 1 and 2, the LO impact factor, 
C'KpI'-.pI'). as in Eq. and 



■'2 J 



A^-^^^{p:,-,k\\k'', 
B^-^^^ip:-,k'(\h 

with z/ = sign(z/a + Ui + 1^2), and 

C^-'^%p+;kt,k-, 

The momentum fractions are defined as 

k 



c^-'^%p:-,k'(\k^,-)A'{h,k2) 
c^-^^\p:-,k^,\k^2')B'ik,,k2) 



(4.3) 
(4.4) 



2C n 



z = 1, 2 [Xi + X2 



(4.5) 



(4.6) 



and the function A'^ as follows: 

with (12) a shorthand for {kik2). Using the dual Ward identity or U(l) decoupling 
equations ||, the function B'^ in Eq. (|4.3|) , and thus the function B^'^^, can be written 

as 

B'{k^, k2) = - [A^(fci, k2) + A^(fc2, fci)] . (4.7) 
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The function C^'^^(p+; fc]*", A;^) is subleading to the required accuracy. The function A'^ 
has a coUinear divergence as 2ki ■ k2 ^ 0, but the divergence cancels out in the function 
B" where gluons 1 and 2 are not adjacent in color ordering 



Using Eq. (|4.7|) , and fixing t ~ — the amplitude ([4 .21) may be rewritten as, 

A^'-''{p:\K\^-\pir,pi') (4.8) 

where the NLO impact factor ioi g* g ^ g g is enclosed in curly brackets, and includes 
six helicity configurations. 

In the multi-Regge limit yi ^ y2, 



hm A^^^np::^ k\\k^2') = C^''{p7; kl^)-C^{q,, k'^^q) 
yi'>y2 ti 



(4.9) 



with qi = —{pa + ki), and ti ~ — and with LO Lipatov vertex, g* g* ^ g [0, ^ 

(4.10) 

Accordingly, the amplitude ( [4.8|) is reduced to an amplitude in multi-Regge kinematics 



C^{q„k\q2) = V2'-^. 

k\ 



47| , with the effective form of a gluon-ladder exchange in the t channel, 



A^'-'''{p':%k'i^\k''2'\piry, 



(4.11) 



with q2 = Pb + Pb' and ^2 



4.2 The NLO impact factor g g* 



qq 



The amplitude g g —>■ qq g for the production of a qq pair in the forward-rapidity region 
of gluon a (Fig. ^) is obtained by taking the amplitudes ( p.ll|) -( p?l3D in the kinematics 



(|4l| ) |4 



A^a-^<^^{p:^,k';\k-''^\pir,pi' 

= 2s 



(4.12) 



Q!2rfi V / d2di 



X ■ 



with ki the antiquark, the NLO impact factor g* g qq in curly brackets, and with 

A^-'^'>{p:-k'^\k^^^) = C^'^\p':^-kl\k:^^^)A'{k,M) 

C^-^\pt;kt,k2) = (4.13) 
C^-'^'^{p:;k^,kt) = 



xlx2 , 
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with momentum fractions as in Eq. ( [4.5D , A'^ in Eq. ( [4 .61 ) and u = Ua- The NLO impact 
factor g* g ^ qq allows for four helicity configurations. 

In the multi-Regge limit kf ^ A;^, the NLO impact factor g* g qq vanishes, since 
quark production along the multi-Regge ladder is suppressed. 



4.3 The NLO impact factor q g' 



qg 



The amplitude qg qg g for the production oi a. qg pair in the forward-rapidity region 
of quark a (Fig. |]b) is obtained by taking the amplitudes (|2.11|) - (|2.13|) in the kinematics 



Pb' ^Pb ) 



(4.14) 



2s U 



J d\a V J d\a 



t 



^9f'''c^''{pl':Pbr) 



with ki the final-state quark, and the NLO impact factor qg* ^ qg m curly brackets. As 
above, the NLO impact factor includes four helicity configurations. 



^ \Pa 1 '^1 ) '^2 y 

B^-^'^%p~^^-k'i\kr^ 

C^'''^{p-;kt,kt] 
C^-^^^{pi-k^,kt] 



C^-'''{p7':k'^X')A'{k,M) 
C^'''{p7'-^l\k',')B'{k,M) 

—iy/Xl 



(4.15) 



1 1 



with A^ in Eq. (^), and given by Eq. (|47|), with v = u-j^. As in Section grT|, the 
function B™ vanishes in the collinear limit. 

In the multi-Regge limit kf ^ /c^ the amplitude ( [4.14| ) reduces to Eq. ( [4.11| ), with 
the substitution ( p.ll|) for the upper line, and the LO impact factor C^'"^ in Eq. (|3^ ). 

The treatment of the amplitude qg — qgg for the production oi a. qg pair in the 
forward-rapidity region of antiquark a is identical to the former, thus the NLO impact 
factor qg* —>■ q g is the same as in Eq. (|4.14|) up to inverting the color flow on the quark 
line |pl|. The corresponding functions A and B are the same as in Eq. ([4.15|). 



4.4 NLO impact factors in the colhnear hmit 



The collinear factorization for a generic amplitude occurs both on the subamplitude and 
on the full amplitude since in Eqs. ( p.l| ), ( |2.11| ) and ( |2.14| ) color orderings where 



the collinear partons are not adjacent do not have a collinear divergence. Hence in the 



^In this context, Eq. (4.7) is only a bookkeeping, since the U(l) decoupling equation is valid only for 
the gluino-gluon subamplitudes corresponding to the quark-gluon subamplitudes used in Eq. (4.14). 
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collinear limit for partons i and j, with fcj = zP and kj = (1 — z)P, a generic amplitude 
( |2.1| ) can be written as 

lim A-'^'--i..., fcf, k"/, ...) = 5^ A-'-{..., P^ ...)Split{7^'^^ , (4.16) 

V 

with / denoting the parton species. Accordingly, for ki = zP and k2 = {1 — z)P, we can 
write the amplitudes (gj), (|]T|) and (|4.14|) as 

= Af'^^^ip:^ 1 ■ spht^f ^^^^(^r, ^2^) , (4.17) 

with as in Eq. ( p.l|) , ( p.4|) and (|3.8| ), respectively, and where we have used helicity 

conservation in the s channel (Section]^). For the collinear factors. Split we obtain 

9^99 (uvi iU2\ 



Split™(A;;'Sfc^^) = tgf"''''^ split's:'' {k';\k. 



2 



>^2 y 



"^2 y 



Split?.-:^''(fc-,fc-) = (7(A=),,,- split!.-:'"^(A;r,A;-) (4.18) 

spiit!.-:'^^(fcr,fc2^) = ^7(A^^),,.-spiiti-;''^(A:r,fc^ 

Splitlr^(A:r,A:r) = ^7(A^^),4split^:-:^"^(A;r,fc^ 
with splitting factors |^ P[ , 

split™(A;+,A;2+) = V2 ■ 

split™(A;i-,A;+) = V2 

split™ (A;+,A;2") = V2 

spliti-^'\kt,k^) = V2 



V4i 


-^)(12) 




z" 




-^)(12) 


(1 


-.)2 


v/^(i 


-^)(12) 


1-z 




(12) 




z 




(12) 





splitr^"'^(A;r,4) = (4-19) 

splitr'^^(fc+,fc2^) = split™(A:+,A:+) = v^^=i-— 

Vl - -2(12) 

spiitr^^(^r,^2^) = spiitr^"^(fcr,fe2^) = V2 

vl — -2(12) 

and split{^'^^'^^(A;]~'^^, fc^'^^) obtained from spXxi^Z^^^^'^ik^^ ^k^^) exchanging (^1^2) with 
[A;2A;i], and multiplying by the coefficient S", Eq. (p.7|), if the splitting factor includes a 
quark pair. 

Summing over the two helicity states of partons 1 and 2, we obtain a two-dimensional 
matrix, whose entries are the Altarelli-Parisi splitting functions at fixed color and helicity 
of the parent parton P] 

^ '^p\Y\.^^^'^\kl\k\^)^p\Yi^p^^'^\k\\kl^)X = 5-'^e'(</'A-0p)p/-/i/2 ^ (4 20) 
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with e^^'^^ '^p^ a phase, where e''^'^+ '^-^ = [21]/(12), and where by definition p^^^^^'^ 
Pl^f'f\ and P^r^'^' = P/+^^^^ and 



'9^99 



P_ 



n9^<l<] 



9~*<iq 



P 



TDl^<i9 



2Ca , 

1 — z 

2Caz{1-z) 

I + (1 - 
z(l-z) 



z 1 — z 

+ + z(l -z) 



(4.21) 



P_ 



1^(19 



9^19 



0. 



l + z^ 



For pi^i9 hehcity conservation on the quark fine sets the off-diagonal elements equal to 
zero, pi^d'i is obtained from pi^'^s Iqj exchanging {z ^ 1 — z). Since we sum over two 
helicity states of the external partons, Eq. ( [4.21j ) is valid in the dimensional reduction 
(DR) scheme |3^. Eq. ( [4.21|) agrees with the corresponding spin-correlated splitting 
functions of Ref. |^ in the DR scheme, after contracting the ones of type P^^f^f^ with 
a parent-gluon polarization as in Appendix ^ The connection of Eq. ( 4.21 ) with other 
regularization schemes (RS) is also given in Ref. |[49|| . 

Averaging over the trace of gq. (jCTp , i.e. over color and helicity of the 

parent parton on the left hand side of Eq. ( [4.20| ), we obtain the unpolarized Altarelli-Parisi 
splitting functions ^ 



2C 



J2 |Split^7^^^^(A;r, 



'^2 y 



S12 



(4.22) 



with C = — 1 for a parent gluon and C = for a parent quark, and where the averaged 
trace of p/^/1/2 is {pf^hf^) = iiPf^f^f^/2 = pIX^'^\ 



5 The Next-to-next-to-leading Impact Factors 

In order to derive the next-to-next-to-leading (NNLO) impact factors, we repeat the 
analysis of Sect. ^ with one more final-state parton. Let four partons be produced with 
momenta fci, ^2, and py in the scattering between two partons of momenta pa and pb, 
with a cluster of three partons, fci, k2 and k^, in the forward-rapidity region of parton Pa, 

yi - y2 ^ y?. '> Vb' ] l^^i^l =^ |A;2±| ^ l^^sxl - K±| . (5.1) 

§Note that in the DR scheme the unpolarized sphtting functions do not coincide with the azimuthally- 
averaged ones. The latter are given in any RS in Ref. Q. 
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Pa a Va h di Vi 




1 ^^lSSlQSLQ^ 

2 voa^mo^l— — 



ny6i5isi5T^mT5i5iy6Ty- ^Tnyw?n5<^n5Tr5T5TP ^w5iy7yi;iW?3Tr5T3^ 



Pb b Vb 



Pb' b' Vh' 



(b) 



id) 



Figure 3: Amplitudes for the production of four partons, with partons ki, k2 and in 
the forward-rapidity region of parton pa- 



5.1 The NNLO impact factor g g* g g g 



We begin with the amphtude for the scattering g g g g g g (Fig. |]a) in the kinematics 
(|5.1| ). Using Eqs. (|2.1|) , (|2.2| ) and (|2.10|) , and the subamphtudes of non-PT type, with 
three gluons of + hehcity and three gluons of — hehcity, |^ , and Appendix 0, we obtain 



s 



(5.2) 

h 



= Ag'^^C^'^{pl^-plY) J2 [A''''{p7;kl?X7.k7^) tr (A'^A'^-iA'^-^A'^-aA^'A 

-X'^\d"yi\d-<y2X'i'"iX^X^' + a^A^'A'^'"3A'^'"2A'^'"i A" - A^'A''A'^'"3A'^'"2A'^'^iA" 
+5<;;39(p^.. k7,\k"a7,kZ') tr [x" X^^^ X^^^ \^ X^^^ 

-A'^A'^'^i A'^'"2A^A^'A'^'"3 + A*A'''A'^''2A'^'"iA''A'^'"» - A^'A^A'^''2_x'^'^i A'^A'^'^a 

with the sum over the permutations of the three gluons 1, 2 and 3, and the LO impact 
factor, C^'^{pI'']p'^Y), as in Eq. ( |3.2|) . From the PT subamphtudes ( |2.1CI| ) we obtain the 
function of ( — h ++) helicities 



M''^ip:-,k'(\k^,\k^,^) = c^-^^{p^^-,k'i\ki\ki^)A'{k,MM) 

where z/ = sign(z/a + 1^1 + 1^2 + 1^3) and 

q± fx^ 1 



A+{ki,k2,k3 



-2 



k 



i± 



X3 (12) (23) 



(5.3) 



(5.4) 



and 



k^ 



1, 2, 3 (xi + X2 + X3 = 1) . 



The functions C^'^^ are a straightforward generalization of the functions C^'^^ defined in 
Eq. (1^ and read. 



^ \Pa ) '>'l 5 '^2 ; '^3 J 



1,2,3 



with z/ = +, (5.6) 
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From the non-PT subamplitudes we obtain the function of ( hehcities 

A9<^9(r,- U+ h~ 



Sl2\k 



Sal2Sl23 
S23 



P{ki, k2, ks)xi H 



X2 X3 



k2, k'ifsai2 7(A:i,fc2,fc3)^gi23|fci_Lp gjaX^Xalg^p 



X2 



S2ZXIX2X3 



S2?,{X2 + 3:3). 



^ \Pa ) '>'l 5 '^2 ? '^3 



Sl2\k 



1± 



■5012-5123 
S23 



, , N , 7(^i,^3,^2)(a;ia;2-a(fci,^3,fc2)(a;2 + a;3))\ f^^. 
-a{ki,k3,k2)xi^ (5.7) 



a{ki, ks, k2)^Sai2 _^ 7(^1, k^, k2)^si23\ki± 
X2 S23a;ia;2X3 

A9'<^9(^- I.- 1.+ ^+^ _ 

^ \Pa ) '^1 ; '^2 ? '^3 

2 



X2 X3 

P _^ Si2xlx2\q±\'^ 

S23{X2 + X3) 



S12 



7(^2,^3, /Cl)^gl23 _ a(fc2,fc3,fcl)'^gal2 ^ ^(fca, ^3, ^1)7(^2, ^3, ^l) -50125123 (a^2 + ^^3) 



5232:1X2X3 



X2IA; 



■3± 



S23X2X3IA; 



3± 



with Sjjfc = (pi + Pj + PkY the three-particle invariant, and 

, ^ _ V^fe3±(v/^gl + v/^[12]) 
a(fci, fc2, ^3) = 

-Sal3 

m,fe2,A;3)^ ^'^" + '^"^[^'J^ (5i 

-Sal2 



„ , , , Va^i X2 X3 [1 2] (1 3) + (2 3) ) 
7{ki,k2,k3) = ^ 



.S123 



as 



Using the U(l) decouphng equations the function B in Eq. ( [5 .21) can be written 

B^''^{p:-k'(\k^,\k^,^)= (5.9) 

IPa ) '^1 5 '^2 ) ^^3 J + ^ IPa ) '^1 5 ^^3 ? '^2 J + ^ IPa ) '^3 ; '^1 ? '^2 JJ • 



In the triple coUinear limit, fci||A;2||fc3, Sect. p.6| , the function A has a double coUinear 
divergence, while the function 5, whose gluon 3 is not color adjacent to gluons 1 and 2, 
has only a single collinear divergence. 

Using Eq. (|5.9|) , we can rewrite Eq. (|5.2|) as 



A^'-'\pi%k\\k^2\k'',^\pir,pl^) = 



(5.10) 
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where the NNLO impact factor g* g ^ g g g is enclosed in curly brackets, and includes 14 
helicity configurations. 



5.2 The NNLO impact factor g g* ^ gqq 

We consider the amplitude for the scattering gg ^ gqq g (Fig. ^d), in the kinematics 
(|5.1| ). Using Eqs. (|2.11|) - (|2.13|) and the subamplitudes of non-PT type, with two gluons 
of + helicity and two gluons of — helicity 13^, we obtain 



1 



= 2sP'^^'^{p:-,k';\k^,\k-''-)- [igf'''C^-^%pl''-y,Y) 
with ^3 the quark, and with NNLO impact factor gg* ^ gqq, 



(5.11) 



\Pa 1 '^l 5 '^2 5 '^3 



Af'^'ip:- kl\k'',\ k^^-) + ( A'^A^'A'^^ ) _ Ar'\p^a: kl\k^,\ k-. 



(5.12) 



+ A'^^A^ A" _ Ar\p7; kVX2\ k^'''^) + A'^^A^A^ _ A^'^P^-, k^k'^^k.. 



d3d2 



d3d2 



+ A'^ A^^ A=' ^ fif ^'^'^ {p:- ; , k^,- ,k^''-)+(X'' X'^ X^ _ B^'' {p^a ; , ^2^ , ^3"''^ ) 

V / d:id2 V / dzd2 

The NNLO impact factor allows for eight helicity configurations. From the PT subam- 
plitudes ( ^.12D -( P?T3D we obtain 

1 



^1 \Pa 1 '^1 ! '^2 5 '^3 J 

^2 yPa I '^1 7 '^2 > '^3 J 

^3 \Pa 1 '^1 5 '^2 5 '^3 J 
^4 \Pa ) '^l ) '^2 ) '^3 J 
\Pa 1 '^1 ; '•'2 ' ^^3 / 

Tjgmqf +. l+ l- 

-°2 \Pa ) '^l ) '^2 ) '^3 / 



-2^^x,a;3-^2) (2 3) 



3^2 ]^ 



fca^ V xi (12) (2 3) 
A:2± 'Va^i (13) (3 2) 



2—^X1X^X3^^3^ (3 2) 
-2 ^ — X2X3- 



(5.13) 



fc3± A;i^(2 3) 

o 9^ 2 1 
-2 - — X, 



A;2±'^'A;i^(2 3) 



with momentum fractions as in Eq. ( |5.5| ). The impact factors from the non-PT subam- 



plitudes p7 



are, 



A9mq(^-.i.+ u- u+\ — 
^1 \Pa I '^1 ) '^2 ? '^3 J — 

2 . ijki, ks, /C2)^gi23 _^ 1/^7(^1, ks, k2)si23a{ki, k^, k2)s3w 



[1 3] (1 2) S23a;iX2X3 S12 (2 3) A;*_LX3yx^ [1 3] k2± 
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ik3± + q±)a{ki, k^fs^w 



\q±\'^xlxy'^ 



^/x^si2\ki±\^k2i_ 



Si2S23|fci±Pa;3 
7(A;i, kz, k2)si23 



X3S23\kl±\'^{x2 + X3) 

Si23y/x^a{ki, ks, k2) 

a{ki, ks, k2){x2 + X3) 



^/X2X3 



-X2 + 



Xl 



(5.14) 



^2 \Pa ) '^1 ) '^2 ) '^3 J ~ 

7(A:i, A;3, A:2)^si23 
[1 3] (1 2) S23a;ia^22^3 



X2X3(-[2 3] V^+[l 3]V^) 



Xl 



S12 [2 3] 



^ k2±{klj_ - [2 3] ^3:2X3)^ ^ VX2X3 ^ ^/X2X3^{ki, k3, A;2)si23 



S23K±Slbb'y/X2X3 



S23 Xl (12) (2 3) [13] A;. 



3± 



(5.15) 



k 3±^/x^oi{k3, ki, k2Ysiw _ 7(^1,^3,^2)^^123 [1 2] 

k2l]k2A]'^S23y/xi X1X2X3 [1 3] S13S23 

(12)a;i Q;(/c3,^i,^2)sifeb'v^((13) - (12)^/x^) 



+ 



(13)s23 (13) |A';2^Ps23v/a^ 

[1 3] ^2X3^^ + (^2^X1 - q]_X2X3)^Jxi 

xr(13) A;*^[2 3] 



(5.16) 



Ar\p--,kt.k2M) = 

{k2± + q±f^/x^o^{kl, k2, k3) {k2± + qA_)qlx\-i{ki, ks, k2)Si23 



(13) \ki^\'^k3^^fx^ (2 3) [13] |/ci±pSi3^XiX3X2 

[12]7(A;i,A;3,A;2)2 S123 , ^1 



[1 3] Si3S23a;ia;2a;3 |A:i_l| 2523(2:2 + 2:3) 
[1 2] a;isi23 [-{k2L + g±)a;2 V'^i^^s [1 3] + S266'7(A;i, A;3, /^2)] 

a;2|/i;i±pSl3S23 [13] 
a:i(fc2± + q±) [1 2] 7(/ci, /C3, fc2)gi23(a^2 + X3) 

Sl3\kll.\'^S23^JXiX3X2 



(5.17) 



Br\p--kt.k2^t) = 

2 ^ -tt(fc3, fci, k2Ysiw ^ {k2± + 5±)^ [1 3] ^/xix3a{kl, k2, ks) 



kl^S23k2±^JX2X3 



\kii_\'^\k3i_\'^k3^y/x^ 



[1 3] ^JXlX3 


;(2 3)v^|A;3x 


' - (fc2± 


+ qi)S23{^\ + 2^2) + ^1±^2± [2 3] X3yx^] 


\f^2 


A;i_L 


^|^3± 


^S23 



|A;i±|2s23V^ 



[1 3] A;i_La;2 + q*x.\/x\X3 



q±X2 
X2 + X3 



+ k 



2± 



(5.18) 
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-(fci± + k2A.){kl^ + [12] ^ h±^2{kl^ - [2 3] v/^)' 



-5235 166' Xg''^ 



X2 



hi. [2 3] S366/X1X2 



\kl±\^\k2±\^S23X3 

+ {kl^ + [1 2] Vil^) ((2 3) gx(fc3\ - [2 3] ./^) - k^^ [2 3] k2±{x2 + x^)) 



\k3±\'^S3bb'XlX2 + \q±\'^\k2± 



2 XiX^ 
X2 + X3 



+ ik*i± + [1 2] y/XiX2){k3^Sibb'X2X3 + fc2±|fc3±P(x2 + 3:3)) | 



(5.19) 



The functions A and B for the remaining hehcity configurations are derived using the 
relations, 



\Pa ) '^1 ) '^2 > '^3 



B 



\Pa 1 '>'l ; '^2 ? '^3 



\Pa I '^1 ) '^3 ' '^2 , 
-°3-i \Pa 1 '•'l ) «'3 5 '^2 



1,2,3,4 
1,2 



(5.20) 



5.3 The NNLO impact factor qg* ^ qg g 



We consider the amphtude q g Q g g g for the production of a quark and two gluons in 
the forward-rapidity region of quark a (Fig. ^) in the kinematics ( |5.1|) . Using Eqs. (p.ll|) - 
( |2.13| ) and the subamphtudes of non-PT type, with two gluons of + helicity and two gluons 
of — helicity [0, we obtain 



A'^'^'''{Pa''\k^i\k'2\k'3' \ p7 ,PI'') 

= 2sl^-^''^\p-''^;k'i\k^^\kl^)\ [igf'''C^-^'{pl'-,p7) 

with ki the final-state quark, and the NNLO impact factor qg* qg g, 

I'-^''\p7';kl\k'',\kl^) 



(5.21) 



(5.22) 



da,, \ d, 



dia 



+ A'^-^ A^-3 Br\p-^^ ■ fcr , k., , ..3 



d\a 



^2 \Pa 1 '^l 1 '^o"2 , fi-o-g ) 



The NNLO impact factor allows for eight helicity configurations. From the PT subam- 
phtudes ( p.l2|) -( pn3D we obtain 

1 



A1'^19g(r)-- U+ U+ 1.+) 
^ \Pa T l^l ^ '^2 ^ '^3 J 



2i 



q± xi 
k 



Br\p^■kt,kt,kt)=2^ 



q± X1X3 



^(12) (2 3) 
1 



k3LS ^2 ki^{2?>) 



(5.23) 



18 



q± xi 



k 



and 



Br'ipl^ fcr, 4, ^3"-) = -^iBf^^'ip^^ kt, kt, kt 
The impact factors from the non-PT subamphtudes 



1,2 



are, 



^ \Pa 1 '•'1 ? '^2 ' ^^3 



2^ / fe2±(feix + [12] y/xTx^)^ ^ 7(/ci, /C2, ksY-fiki, ks, k2fs 



2„2 
123 



|A;i±ni2]s36b'y^ ' (xiX2X3)3/2si2S23 (2 1) [1 3] 
[1 3] (1 2) Xi(^/X^(2 3) - + ^3, /C2)S123 



S23|/^i±Pa;3 



Sl2S23^1±A/a^ 



2 2:1X2X3 



X2 + X3 



+xi(-g±/c3_LX2 + A;2±a/x3(^i±v^ - [2 3] V^)) } 



^ iFa ) '^1 ) '^2 ) '^3 J ~ 

2i 



(3{ki,k2,k^fs^wk*2^ ^ Pjki, k2, k'i)sm'k*2^y/x{{xi. + X2) 



[12]si2|fcixPxf 



si2\ki±\'^y/x-iX2 [3 2] 
fc;_Lxf^ [g^7(/ci, A;2,/i;3)*Si32 



Si2|fcixn3 2] 



+ Si32a/x^ 



%,l3{ki, k2, k3)ssbb'ql 



Xiy/xsk 



2± 



2„2 



+ 



'132 



7(fcl,fc3,fc2)*7(fci,fc2,fc3; 
(xiX3X2)3/2si2S32 (3 1) [1 2] 



(5.24) 



(5.25) 



^ 3^1 [-gi2|g±Pa^2Va^i3^3 + (13)/3(A:i, A:2,fc3)g3feb'fe2±(^3 + 3^2)^] 
Si2S32|/i;i±Pa;2v/^(a:3 + ^2) 

[1 2] (1 3) k;^^/4F2 + P{ku k2, k3){k*2^Si32^/^3 - (1 3) kl^qlx2) 



[12] Si2S32|^l±IV^2X3 



(5.26) 



Br^{p-;kt,k2,kt) = 

^. f {k2± + q±y [l^]xiy/x^ ^ A;2±[13]xiX3 a{k2,ki,k^)*a{k:i,ki,k2fs\^^, 



+ 



|/i;i±|^^3±S266' 

[13]y^ 



|^l±|^|fc3±|^S23X2 
3^3 



|A;i±|2|/C3±|^S23X2v/^ 

* / qi\k2,i_\'^x\x2 



(2 3 ) A;*^ I fc3± 1 2 ^ I A;3x 1 2 A;*^ /ca^ S23 v^^xs 

[-|^3±|^(fc2± + g±)a;iX2 + |A;i±|^A;2±(a;i + X2)x3] 
k2Lkli_\k^i_\^xi + fc2±|fc3±|^ [1 2] x^l'^^Jxi 



X2 + X3 



+ |A;i±| A;2x(Xi +X2)X3 



(5.27) 
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T^qigggf -. l+ l+ u- 



2i 



a{k3, ki, k2)*a{k2, fci, ksfslf^,,, (1 3) y/x^a{ki, k2, ksfsibb' 



k3±k2±\k3±\'^S32^/x^X2 



\ki±\'\k3±\'k 



3± 



|^l±n^3±P-523 



[1 2] (1 3) — ^^(s2feb'Xi - 3^1^x2) 



+ V^ii-ki±q±xr + (1 3) k3±y^iix3 + X2)) { qf - [1 2f - ^^-^gl [1 2] 



3/2 

-r 



X3 + X2 



+ (13) {Siw + S2bb')^/xi 



Br^ip--kt,k^,kt) = 

2i r {k2± + q±)xi{kl^ + [1 2] fc2±(fci\ + [1 2] v^ili^)' 



l±r L ^3± [12] ^/x^ 

{k2± + q±f[l3] xi^ 



[12] S3bb'y/X^ 



k3±S2bb' 

Br'^ip-;kt,kt,k^) = 

2i J (1 3) y^a{ki, k2, k3fs2w_ ^ Piki, k2, k^fs^bb'k^^ 



S12 [2 1] X2 



3/2 



Xl 



{12)x3a{ki,k2,k3)s2bb' 



Xik, 



3± 



(12) k;^ 

, A , gl^i - k*^X2 \ Piki,k2,k3)s3bb' 1 1 
V V^[12] J ^2 J J 



(5.28) 



(5.29) 



(5.30) 



5.4 The NNLO impact factor qg" qQQ 

We consider the amplitude q g ^ qQ Q g for the production of three quarks in the forward- 
rapidity region of quark a (Fig. ^d) in the kinematics ( |5.1j ). Using Eq. ( |2.14|) and the 
subamphtudes of non-PT type, with two gluons of opposite hehcities [Q, we obtain 

A'''-"^^^'{p-''\k'^\k^,\k^''- 1 = 

1 



2s P-"^^ {p-^^ ■ k\^ , k!;- , k^-"- ) - tg f'^C^-^^ ipl" ; p^r ) 



with NNLO impact factor qg* ^ qQQ 
P''>^Q{p-''^;kl\k^,\k^^-) = 



(5.31) 
(5.32) 



^dsa^dida^l {Pa ) 5 '^2 ) % ) Tf^dia", 



• •aia^d3d2 2 



a ) ft-l ) "-2 > "^3 i 



Nr. 



-5gQ (1^3). 
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The term proportional to 6qQ is due to the interference of identical quarks (i.e. with the 
same flavour and helicity ) in the final state. The NNLO impact factor allows for four 
helicity configurations. From the non-PT subamplitudes we obtain, 



Ai-^qQQ(„~. u+ u- u+ 



+ 



X2 





(1± 


2 


1 - Xi S23 




2 



^/x^S23\kl±\ 



+ 



Xiaih, ki, k2) Sa23 - 



+ 



l{ki, ks, k2 



7(fcl, fcg, fc2)Si23 



^2 \Pa 1 '•'I ) '^2 5 ^^3 J 



V^2a^S23 



7(fci,fc3, A;2) 



+ ^/x^a{k3, ki, ^2) 



3^3 
Xi 



^/xlx2X2. 1 




2 


l-Xi S23 


\ki± 


2 



+ 



Xi a{k3,ki,k2 



a/^2^"S23|A;i± 



\k 



l-L 




Xia(A;3, ki, ^2) 5^23 



a;2a;3 
7(A;i, A;3, ^2)5123 



■§23 



(5.33) 



(3{ki, ^3, ^2) + ». —a{ki, k^, /cg) 
3:3 V 2^2 



with a, /3, 7 defined in Eq. (|5.8| ) and 

/ig;'?OQ/ +. z^- ^- \ _ A<j;qQQ ( - . 1.+ t.- 1.+ 

\Pa 1 '^1 5 ^^2 5 ^^3 J ~ \Pa i '^1 ) «'3 ' ^^2 



j^Q,qQQ , k2 , k^ 



b: 



Q;qQQ 



i = l,2 
(Pa 'i kf 1 k^ , A;^) i = 1,2 . 



(5.34) 



Note that for each helicity configuration, we have the following relation between the 
functions A and B 

(5.35) 



j^<r,qQQ _j_ j^mQQ _ j^^mQQ _|_ 



5.5 NNLO impact factors in the high-energy Umit 



The amplitudes (|5.10| ), ( ^.11|) , (|5.21|) and ( |5.31|) have been computed in the kinematic 



limit (|5.1|) , in which they factorize into an effective amplitude with a ladder structure, 
made of a three-parton forward cluster and a LO impact factor connected by a gluon 
exchanged in the crossed channel (Fig. In the limits yx — y2^ Vs or 2/1 S> ?/2 — 2/3, the 
amplitudes must factorize further into NLO impact factors or into NLO Lipatov vertices 
for the production of two partons along the ladder. Such limits constitute then necessary 
consistency checks, and we display them in this section. 

In the limit, yi ^ 1/2 — ys, the NNLO impact factor, g* g ^ g g g, Eq. ( p.lO| ), factorizes 



into a NLO Lipatov vertex for the production of two gluons convoluted with a multi-Regge 
ladder (Fig. |a) 

hm I {tgf Yl r'-^''r'-^'''r''-^''"A^-''^{p:-, c% O) 1 (5.36) 
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Figure 4: Limits of the amplitude for the production of three gluons in the forward- 
rapidity region of gluon pa, for yi ^ ^2 — Us (a) and yi ~ ^2 ^ Us (b). 



with the NLO Lipatov vertex, g* g* —^g g, for the production of two gluons k2 and 
given by [0 g 



k2± [ S23 
_^ {q2± + k3i_Y _ q2± + k3± [/c^ + k~^ 



P ^ kl^\q2±\ 



(5.37) 

S3bb' k2±k3j_ 



(^2 "I" ^3 )^3 (^2 "I" ^3 



^2 ^3 



■536b' 



S23 



k2± - 



^2 "t" ^3^ 



'3± 



with exchanged momenta in the t channel qi = —{pa' +Pa), 92 = Vv + Ph, three-particle 
invariant 536^ = (^3 + ^2)^ — — (|'?2± + ^3±P + ^2^^^); and with the mass-shell conditions 



k: 



\hi?/kt for i = 2,3. 



In the coUinear limit, /c2 = zP and /cs = (1 — z)P, the NLO Lipatov vertex ( ^.37[) reduces 
to the splitting factor ( [4.19| ), and amplitude (|5.10|) factorizes into a multi-Regge amplitude 
(|4.11| ) times a coUinear factor ( [4.1^ ) 

k2\\k3 



Y.^''^''(Pa^ ^1 I P'^ I p7^P?) ■ Splitl7''(fc2% ^: 



V3\ 
3 )i 
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In the limit, ?/i ~ ?/2 ^ I/3, the NNLO impact factor in Eq. ( p.lO| ) factorizes into a NLO 
impact factor, g* g ^ gg, Eq. ([4 .81) , convoluted with a multi-Regge ladder (Fig. Ijb) 



0-G53 



(5.38) 



0"2 . 



1 



with gi = — (pa + ^1 + ^2)5 and with LO Lipatov vertex C^(gi, /cg^, 52), Eq. ( |4.10| ). 

In the limit, |/i 3> 1/2 — Vs, the functions A and B in Eq. ( |5.13D -( |ODD fulfill the 
relations Af^^"" = A^'^"^^ = 0, Bl''^^" = -Af'^", and Ef'^^"'? = thus the NNLO 

impact factor, g* g ^ gqq, Eq. (|5.12| ), factorizes into a NLO Lipatov vertex for the 
production of a qq pair convoluted with a multi-Regge ladder (Fig. ^), 



lim F^^^"^(p,^^A:r^A;2^^fc3"' 
j/i>j/2-y3 



(5.39) 



(A^'A^),3,-,A^"^(gi, A;3-^^ g2) + K"'. k^^q,) 



with the NLO Lipatov vertex, g* g* ^ qq, for the production of a qq pair ||23|, pi 



-2 




+ 



^3 ^3±ki±l 



^3 1(^2 +^3)'523 k2±ik2 + ^^3 )S23 



+ 



^3 ^2±(?2± + fc; 



3±J 



^2 ■^366' 
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with gi, q2, and Ssw as in Eq. ( |5.37| ). 

In the colhnear hmit, k2 = zP and = {1 — z)P, the NLO Lipatov vertex ( ^.401) reduces 
to the sphtting factor ( [4.19| ), and amphtude (|5.11|) factorizes into a muhi-Regge amphtude 
(|4.11| ) times a colhnear factor ( [4.1^ ) 



= k\^ I I ■ spht^-«"''(fc2^^ ^3^^) . 



In the limit, 1/2 — Us ^ Vi, the functions A and B in Eq. ( p.l3p -( pr^ ) fulfill the relations 



Af 



am 



0, 5f 



999 



-Af^''^ and Bl 



ami 



-Af'^''' thus the NNLO impact factor, 



Eq. (|5.12| ), factorizes into a NLO impact factor, g* g —>■ qq, Eq. (|4.12|) , convoluted with a 
multi-Regge ladder (Fig. |b). 



hm F'^^'^ip:-,k'(\k!^\k^' 
j/2=;j/3>s/i 



with qi 



ti 

-{Pa + k2 



^gr''"''C^{qiX\q2) 

h). 



(5.41) 
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In the limit, ?/i 3> 1/2 — Uzi the functions A and B in ( |5.23| )-( |5.3CI| ) fulfill the relations 



f>q;qgg (^-1^1. u^i u>^2 u'^3\ _ j^gmg (^-1^1. ui^i 1,^3 i.^2\ 



thus the NNLO impact factor, qg* Qgg, Eq- (|5.22|) , factorizes into a NLO Lipatov 
vertex for the production of two gluons (|5.37|) convoluted with a multi-Regge ladder 
(Fig. |a) 



US'- 



(5.42) 



■0-3 ) 92 J 



In the limit. Hi — 1/2 ^ Us, the functions A and -B in ( |5.23| )-( [5730| ) fulfill the relations 



^ \Pa I '^1 5 '^3 5 



l'2\ 
2 J 

T^q;qgg ( ^-yi . h.yi 1.1^2 h.v3\ 
T^q;qgg ( ^-I'l . uvi lus iu2\ 

^2 \Pa 1 '•'1 ? '^3 5 '^2 J 



_ Aq<qgg(r,-i^i. u'^i 1.^2 h.T^z\ 

^ \Pa 1 '•'1 ? '^2 ? '^3 J 

f?q;qgg(r,-'^i- u"^ h"^ h''^\ 
\Pa 1 '^i 5 '^3 5 '^2 J ' 



thus the NNLO impact factor, Eq. ( 5.22 ), factorizes into a NLO impact factor, qg* ^ qg, 
Eq. ( [4.14|) , convoluted with a multi-Regge ladder (Fig. ||b) 



hm /'?^^^^(p-'^-fcr^A;2^^A;3^^) 



(5.43) 



9 



tgr''^'' C^{qukl\q,) 

ti L 

with qi = -{pa + ki + /cs). 

In the limit, yi ^ y2 — 2/3, the function A2 in Eq. ( |5.33|) vanishes, A^''^^'^ = 0, and 
using Eqs. (|5.33| )- (|5.35| ) the NNLO impact factor, qg* qQQ, Eq. (|5.32|) , factorizes into 
a NLO Lipatov vertex for the production of a qq pair (|5.40|) convoluted with a multi-Regge 
ladder (Fig. 0) 



A'-'^'iPa"'; fcr, ^2^) + (a^'a'^^) B^'''^ip-''^- k\\k 

i\a \ / dia 



hm p-''^Q<^ip^''^;k^,\k^2',K'"'] 

j/i>y2=i|/3 



[9X,,-aC^''ip-''';k^,' 



(5.44) 



(A^'A^),3,-,/l^"^(gi, A;-, ^3'^, g^) + {XX^^.^.A^^iqu A;-, ^2) 
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Pa a Ua h dj, 




Pb b Vb Ph' b' Vy 



Figure 7: Limit of the amplitude for the production of a quark and a quark-antiquark 
pair in the forward-rapidity region of quark p^, for ^ ^2 — Z/s- 



5.6 NNLO impact factors in the triple coUinear limit 

In the triple collinear limit, ki = ZiP, with zi + Z2 + Z3 = 1 a. generic amplitude must 
factorize as 



11 



lim ^■••'^i^2'^3. 

fclp2||fc3 



^2^) ^3^; ■••) ~ ^ 



(5.45) 

Accordingly, we must show that taking the triple collinear limit of the NNLO impact 
factors, we can write the amplitudes ( |5.1(]| ), ( |5.11| ), ( ^.21| ) and (|5.31|) as 



A:il|fc2p3 



Pw'^^^Ph) 



^3 J 



and (13^ 



with / denoting the parton species, J[-^9^f9 given in Eqs. (|3.1| 
Split-^^-^^-^^-^^ the polarized double-splitting functions. 



(5.46) 
and with 



In the triple collinear limit, the functions A of Sect. f).l[ |5.2|, |5.3| and yield a 



quadratic divergence as S123 ^ or Sij — > with i,j = 1,2,3. In the same limit, the 
functions B have a single collinear divergence since only two out of the three partons are 
color adjacent. However, terms with a single divergence when integrated over the triple 
collinear region of phase space yield a negligible contribution |]lOl, thus we ignore them. 

It is easy to show that a function A^'^^, Eqs. ( ^.3|) -( pl7|) , differs from its reflection by 
a term which contains only a single divergence. Using this property and Eq. ( |5.9| ), we 
obtain a dual Ward identity and a reflection identity for the functions A^'^^, up to singly 
divergent terms. 



^ \Pa ) '^1 5 '^2 ) '^3 J + ^ \Pa i '^1 ) '^3 ? '^2 J + ^ \Pa ) '^3 ' '^1 ' '^2 J ~ ^ 

A'^'^'ip:-, k'i^k'i^kl^) = M''%p:-, kl\ fcD . (5.4?) 



Using the identities ( |5.47| ) in the impact factor in Eq. ( ^.10| ), we can factorize the color 
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structure on a leg 



CT3 N 



where (F'^)^^ = ^/'"'^ Thus amphtude ( |530D can be put in the form of Eq. ( [5.46| ) with 
coUinear factor 

spht!.-:''(A;r^fc2'^fc3^^) = (F'^-F'^-)^,^spht?.-:''(CsA;:-,fc3^«). (5.49) 

The sphtting factors spht^!.'^'^^ are the functions A, Eqs. ( |5.ij| )-( p?^ ), in the triple collinear 
limit, up to singly divergent terms, and thus they fulfill the identities, Eq. ( 5.47| ). The 
splitting factors of PT type can be soon read off from Eqs. (|57^ )- (|5.6| ), while for the ones 
of non-PT type we note that the coefficients of Eq. ( ^.8|) reduce to 



zi + 2:3 



(5.50) 



with 



Thus we obtain 



7(fci,A;„fc3)- ^^^5(1,2,3) 
■S123 



5(1, 2, 3) = [1 2] (v/iT (1 3) + v^i^ (2 3)) . (5.51) 



v^(12)(23) 



splitr'^(A;+,A:2+,A:^) = 2—^ , }, , 
splitr''(A;r,fc2^,fc3+) = 

splitr^^(A;+,A:2+,A:o-) = 2—^- — - — - (5.52) 



2 



..... + 3)^(1,2,3) 



splitr'^(A;r,A:2+,A:3+) = split?,-^='^(A;3+, ^^2+, A^f) 

spiitr"^^(fc+, /ts", k+) = -spiitr'^^(fc2-, kf, k+) - spiitr"^''(fc+, k+, k^ 



27 



In the triple collinear limit of the NNLO impact factor g g* — > qqg, the functions 
A^a^'J, Eqs. (p7[3|)-(|5:20|), fulfill the relations Af'^^"^ = -Af'^^"'^ and Af''^'^" = and 
^f^'^'^ba"; K\k'2',k^'"') = Af^^'^p'^^'^; kl\k^''\k''^^). Thus amplitude can be put in 

the form of Eq. ( |5.46|) with collinear factor 

Split™(A;r^A;2^^fc3-^^) (5.53) 



jA'^A^^),^,- split™(A:r, A:-, ^3"-) + (A'^^A^),^,-^ split™ (fc-, ^3"-, k 
with 

1 



split™(fc+,A:2-,fc3'") = 2i 



^1 (12) (2 3) 



spht™(..^4,^3)-2.3y|^^ (5.54) 



split™(A;+,A;2-,fc3+) 



■5l2'523 



• ^(1,3, 2)^ [12] 
[1 3] S123 



(5(1, 3, 2) ^ 22(-2l - 2:3) [2 3] ^^ ^ v^(-22Sl3 + ;Z3S23 + ZiZ2Si2'i) 



z?, V \f^\ [13]23 / v^(l-^i) 

(5(1, 2, 3)^ [13] (5(1,2,3)(1-Z3) V^gi2 
[12]si23 (l-^i) 



split™(fc+,A;2+,fc3"): ^ 



S12S23 



Writing the functions A, Eqs. (|5.23|) - (|5.3CI| ), in the triple collinear limit of the NNLO 
impact factor qg* ^ qg g-, the amplitude ( |5.21| ) can be put in the form of Eq. ( ^.46|) with 
collinear factor 



Split™(A;r, fc2^^ kf) =g'J2 (A''-A'^-3)^^^_split™(A;r, C^, C^) . (5.55) 

with 



5;(12)(23) 



split™(A;+ A:2+,A:3+) = - 
spht™(.r,4,^3^)-^^^ (5.56) 



split™ {k+,k^,k 
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2i 



S12S23 



5(1, 3, 2)2(^ [1 3] + [2 3]) ^ 7^^(1-^3)^(1,3,2) ^ 711^2^12 +^^33) [1 3]' 



[1 3] S123 1-2:1 
2? 

splitr''^^(A;+, A;2+, A;3-) = x 

■Sl2'S23 

5(1, 2, 3)^(7^1 [1 2] -y/j^ [2 3]) ^ 7^^(1-^3)5(1,2,3) _^ v/irz2g23 ' 

[1 2] S123 1 - ^1 , 

In the triple collinear limit of the NNLO impact factor qg* —>■ qQQ, the functions A 
dPI) fulfill the relation Af"^^^ = Af^^^ . Thus the amplitude (lOTD can be put in the 
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form of Eq. ( p.4(j| ) , with collinear factor 



Split 



q^qQQ 



(h^ h^^ h^'^ 



(5.57) 



where the second term occurs for the case of identical quarks, and c is the color index of 
the parent quark. The splitting factors are 



split 
split 



S23 

q^qQQfu- u- u+ 



y/ZlZ2Z3 _^ (5(1,3,2) 



' (^1 ? ^2 ? ^3 



1 - Si23 
Splitr^««(fc+,fc3-,fc+). 



(5.58) 



The factor sp\itl~^^^^^^''^{ki'^^ , k2'^^ , k^'^^) can be obtained from split 



in Eqs. ( |5.52| ), ( p.54| ), ( |5.56D and ( [5. 58] ) by exchanging (ij) with [ji], and multiplying by 
a coefficient S, Eq. (p.7|), for each quark pair the splitting factor includes. 

Using Eq. (^M) and Eqs. (|532D -(F58D, and summing over the two helicity states 
of partons 1, 2 and 3, we obtain, as in Section the two-dimensional Altarelli-Parisi 
polarization matrix at fixed color and helicity of the parent parton. 



^ Split{^^^^^^^(A;;'\ ^2^% fc^^)[Split^^^^^^^^(A;r, '^r, '^D]* = 5'"' ^P^;;^'^'^' , (5.59) 



'123 



where P_ 



P 



and P 



. For splitting functions of type 



pq^qf2h^ namely for pi^i99^ pq^qQQ ^nd P'/i^9i92g2^ ^^ere the last splitting function is 
for identical quarks, helicity conservation on the quark line sets the off-diagonal elements 
equal to zero. 



Averaging over the trace of matrix ( |5.59| ), i.e. over color and helicity of the parent 
parton, we obtain the unpolarized Altarelli-Parisi splitting functions [|10| 



^ J2 |Split{7^^^^^^(A;-,fc-,^, 



4g' 



(Pf- 



*fihh 



(5.60) 



'123 



with C defined below Eq. ( [4.22|) . For l^pa^a^d^g-i^ ^ ^]-^g ^\xm over colors can be immediately 
done using Eq. (|2.5| ), and it yields 



|gplits_-gi^;2g3(^.i^^.2^^.3)|2 ^ 4c^(^j ^ \sp\ii'_-:^\kZ\kZ\k 



I'.raN ,2 
0-3^1 • 



(5.61) 



<tG53 



with Cn{Nc) as in Eq. ( p.8|) . Eq. (|5.61|) shows that for the purely gluonic unpolarized 
splitting function the color factorizes. 

Since the averaged trace of jg (^pf-hf2h^ = ti pf^f^f^f^ /2 = pI^^^^^^\ we 

have checked that for the diagonal elements, P^'^^^^'^^''^ , our expressions agree with the 
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unpolarized splitting functions of Ref. |]TU| by setting there the RS parameter e = 0. 
Finally, for the off-diagonal elements of the splitting functions of type P^^s-^^/a obtain 



P: 



■S123 



^(TlO-3 



■^O"! 0-2(1 ^a'3)^(T3 



2(1 - z„2)z, 



(1 20-3)^0-3 



(1 — ^0-3) ^0-3 

^(71 ^0-3 -Dg-g ( 1 2zo-2^ 

(1 ^0-2)^0-2 



with 
and 



Di = [i j] y/zj + [i k] y/zk with i,j,k = 1,2,3 and j,k i , 



Pi: 



- ( (7^P^^5'i'^2g3 (ab) _|_ ^ ^pg^giq2q3 (nab) 



where the abelian and non-abelian terms are, 



P: 



9^919293 (ab) 



2s 



123 



+- 



5l2Sl3 

19^919293 (nab) \ ^ ■^123 



E 

f7G52 



1 

~2 



^1-^1 + V 2 ^f3 -0(T2 -0(73 "I" 2 



'1 - 2;i)2;i 



2D2 J2 ^ 

_ n2 _1_ '^3^0-2^°-3 



(5.62) 



(5.63) 
(5.64) 



(5.65) 



We have checked that Eqs. (^1621)- (^1651) agree with the corresponding spin-correlated 
splitting functions of Ref. ||lT| after contracting them with a parent-gluon polarization as 
in Appendix 0, and after setting the RS parameter e = 0. 



6 Four-Parton Forward Clusters 

The procedure of Sects. ^ and | can be clearly extended to n-parton forward clusters. In 
a forward cluster there are one incoming and n outgoing partons. Thus, for purely gluonic 
clusters there are 2"'"'"^ helicity configurations. However, in the high-energy limit two of 
these are sub leading, thus an n-gluon forward cluster contains 2(2" — 1) helicity configu- 
rations. For n-parton forward clusters including qq pairs, all the helicity configurations 
are leading; then an easy counting yields 2" helicity configurations for the one including 
a. qq pair, 2"~^ for the one including two q q pairs, and so on. For n = 3, we obtain the 
helicity configurations dealt with in Section ^. 
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(d) 



Figure 8: Amplitude for the production of five gluons, with gluons ki, k2, k^ and k^^ in 
the forward-rapidity region of gluon pa- 



6.1 The NNNLO impact factor g g* g g g g 

Here we analyse in detail the four-gluon forward cluster. We take the production of five 
gluons with momenta /ci, A;2, /ca, k^ and py in the scattering between two partons of 
momenta Pa and p^, and we take partons fci, A;2, k^ and k^ in the forward-rapidity region 
of parton pa (Fig. ||a), 

1/1 - 1/2 ^ 1/3 - 1/4 > w ; \kiA-\hi-\-\hA-\kii-\-\Pb'A- (6-1) 



Using Eqs. (^^ ) and (|2.1CI|) and the subamplitudes of non-PT type, with four gluons 



of -|- helicity and three gluons of — helicity we obtain 



a''^''{p:% ^rs k^K' i k^k^) = (6-2) 
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tr ( A"A'^-iA'^'^2A'^-3A'^-4A'''a* - A'^A'^^iA'^^2A'^'^3A'^-4A^A^' 
-A^A^'A'^"4A'^'^3A'^'^2A'^'^i A" + A^'A''A'^"4A'^'^3A'^'^2A'^'^iA''j 



+ 



ya. hy<^l ]/'<^2 h.^'^S i''^'^4^ 



tr ( A"A'^'^iA'^'^2A'^'^3A'''a''A'^^4 _ a"A'^^iA'^^2A'^-3A''A'''A'^^4 



-A''A'''A'^"3A'^'^2A'^''i A^A"*"* + a^'A^A'^"3A'^'^2A'^'^iA"A'^"4 j 
ind'^^afn'^a. u'^"! hy<^2 h.^'^i ty'^A\ 



with the sum over the permutations of the four gluons 1, 2, 3 and 4. From the PT 
subamphtudes ( p.lO| ) we obtain the functions of ( — h + + +) hehcities 



^ k2^, k^^, k^^'^) — C^'^^(Pa"! ^1^5 ^2^? ^3^' k4^'^)A'^(ki, k2, k^, k^) , 
where z/ = sign(z/a + ui + 1^2 + t^s + 1^4) and 

1 



(6.3) 



and 



A^{ki,k2,k3,ki) 
kt 



As in Eq. (p.6|), the functions C^'^^s are 



ki±y X4 (12) (23) (34)^ 



1,2,3,4 (xi + X2 + X3 + X4 = 1) . 



(6.4) 



r^'^^d ( r^a . h.Vl h.V2 I.U3 T 



= - 



X2 Ui 



I = 1,2,3,4 



with u = +, (6.5) 



From the non-PT subamphtudes |^ we have obtained the functions of ( h + +) 

hehcities. We do not reproduce them here because they are quite lengthy. They are 
available from the authors upon request. 

Using the U(l) decoupling equations for one and two photons, the functions B and D 
in Eq. (|6.2| ) can be written as 



"^4 y 



^ \Pa 1 '^1 5 '^2 5 '^3 ) '^4 

\Pa ) '^1 5 '^2 ? '^3 ' '^4 / ^ \Pa ) '^1 ) '^2 ? '^4 5 '^3 / 
^ ^ \Pa 1 '^1 ; '^4 ? '^2 ? '^3 / ^ \Pa i '^4 ? '^1 ' '^2 ' '*^3 )\ 



(6.6) 



"^4 , 



"^4 y 



ngAg ( ^Va . I.U2 I.U3 ih 

^ \Pa 1 '^l 1 '^2 5 '^3 ? '^4 

yPa 1 '^1 ) '^2 ? '^3 ' '^4 J ^ iFa i '^1 ? '^3 ' '^2 ' '^4 

^ ^ \Pa ) '^3 5 '^1 ) '^2 ) '^4 J ^ \Pa i '^1 ) '^3 ' '^4 i '^2 I 

_|_ AgAg ( . 1.^3 UI^l 1.1^4 U'^2\ < A g-Ag(^l^a. 1.1^3 1.1^4 l.l'l l.'^2\] 
^ ^ \Pa 5 '^3 5 '^1 5 '^4 5 '^2 ) ^ ^ \Pa ) '^3 ' '^4 ) '^1 ? '^2 JJ • 



(6.7) 
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In the quadruple coUinear limit, A;i| |/c2| j/cal |/c4, Sect, p^ , the function A has a triple 
collinear divergence; the function B, whose gluon 4 is not color adjacent to gluons 1, 2 
and 3, has only a double collinear divergence; the function D, where gluon 1 is adjacent 
to 2 and gluon 3 is adjacent to 4 but the pairs are not adjacent one to another, has two 
single collinear divergences. 

Using Eqs. ( |6.6|) and ( p.7|) , we can rewrite Eq. (|6.2| ) as 



A9 9^5g( Ua u^i 1.^3 I ^^b' y 

\Pa 5 '^1 5 '^2 ? '^3 ' '^4 I PV 1 Pb 



(6.8) 



2s { {igf r'^'^i7^'^-2'=7'''^''^'7''''^"*''''^^'^^(Pa^C\C%fc-?,fca 



where the NNLO impact factor g* g gggg is enclosed in curly brackets, and includes 
30 helicity configurations, in agreement with the counting above. 



6.2 NNNLO impact factors in the high-energy hmit 

The amplitude ( |6.8| ) has been computed in the kinematic limit ( |6.1| ), in which it factorizes 
into a four-gluon cluster and a LO impact factor connected by a gluon exchanged in the 
cross channel. In the limits ~ ?/2 — I/3 ^ I/4 or ?/i ~ ?/2 ^ I/3 — I/4, or ?/i ^ ?/2 — I/3 — I/4, 
Eq. ( |6.8|) must factorize further into a NNLO impact factor or into a NLO impact factor 
times a NLO Lipatov vertex, or into a NNLO Lipatov vertex (Fig. §), respectively. While 
the first two limits constitute necessary consistency checks, the last one allows us to derive 
the so far unknown NNLO Lipatov vertex for the production of three gluons along the 
ladder. 



In the limit, ?/i ^ ?/2 — I/3 ^ Ua, the NNNLO impact factor, g* g ^ g g g g, in Eq. ( |0| 
factorizes into a NNLO impact factor, g* g ggg, Eq. ( p.lO| ), convoluted with a multi 
Regge ladder (Fig. 



3' 



1 

X — 



zgr"'-^"'C\q,,kl\q2) , (6.9) 



with qi = —{Pa + ki + k2 + k^), q2 = Pt' + Pb, and with LO Lipatov vertex C^(gi, ^3^, ^2)5 
Eq. ( CT . 



In the limit ?/i ^ ?/2 ^ I/3 — Va, the NNNLO impact factor in Eq. ( |6.8|) factorizes into a 
NLO impact factor, g* g g g, Eq. (|4.8|) , times a NLO Lipatov vertex for production of 
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yi-2/2 



two gluons g* g* g g ( |5.37D , convoluted with a multi-Regge ladder (Fig. |]b) 

^ I iwY E /^''^-3^7^"'^.4^"'A«.(g^, C% A;:?, 92) ^ , (6.10) 



with qi = -{pa + ki + /C2). 



6.3 The NNLO Lipatov vertex 

In the limit ^ ?/2 — 2/3 — 2/4? the NNNLO impact factor in Eq. ( |6.8|) factorizes into a 
NNLO Lipatov vertex convoluted with a multi-Regge ladder (Fig. ||c) 



crG04 



[igr'^'C^-^'ipZ-.K 



(6.11) 



t 



q2) 



with the NNLO Lipatov vertex, g* g* — > ggg, for the production of three gluons k2, k^ 
and ^4 enclosed in curly brackets in the right hand side, with 



A'9^q,,kt,kt,kt,q2) = -2V2 



X2 



(ii± q2± 



X4 (2 3) (3 4) k2± 



(6.12) 





\(ll±\ 


\'k2AKJHq;^ +K± 


)X2 


S4bb'k2j_ S23 {\k3±\ 


P X2 + 







k 



2± 



M± Q*2± K± (3 4) + (2 4) (g*^ - k;^ ) {q;^ + k*,^ ) [2 3]) ^ 



S466'(3 4)A;^_L S23 

-Kl HqIl - ) S23 + q*ii_ S34 ) 







P /^2±(A;3^)^a;2 


' (34)/c*^ , 


§23 




p 2:2 + 


|/^2± 


P X3 ) ^/xT 



-q*i± {-k2±k*^ 7x7+ (2 3)A;|^ + ^2±[2 3]v^ + (2 4)A;*^ 



34 



+ k*,^ {{{23)kl^ +k2±[23])^+{24)k*^^))) 
^ \q2±\^X2 ((2 3)Vxr+(2 4)V^) ^ g2±_fc2±(?i\ - k*^ f {q*^ [34]^) 



(2 3)(34)s234^/^ 

(M^^3± kl^ 



S2ao'S466'(34)A;2± 



(34)s23 i\qi±-q2±\ + ^234) (|A;3±| + \k2±r ) ^ 

\qi±\^ fe2x[3 4] ((2 4)fc:^ y^+{23)kl^ y/xT) 
X3 X4 (34)s23 S234 (|gi± - g2±P + S234) 

[34] 



qt± k2±X3 X4 ((2 3)Vir+ (2 4)^^) 



(34)S23 S234V^^3 X4 

+q2± Vx^ik4± {q2± + kl^ ) X2 x^ + (|A;3±P a;2 + k2±ik^^ - ql^ )x3 ) X4 )) } (6.13) 



A^^{q,,k+,k^,k+,q2) = 2V2 X 

" |g2±P fe3± v^v^ _^ |g2±P (gl± - fc2±) A:^_L 3^3 _ |g2±P \/3:2 X^ 
_(2 3)/C2±S34 (1 - a;2 ) S2aa'^2±S34 (1 - a;2 ) (2 3)S34 {I-X2) 





?1± 


' (^3± )'fc2\ 


(92± +kl^ 


)X2 


S4bb'k2±S23 ( 


/i;3± 


2 X2 + |/i;2± 


'xs) 



{k3± fV^ (- {q*2i. + A;:^ ) [2 3)^x7 + gi\ (g2\ + [2 4] V^)) 

S466'^2±S23 

q*i± k3± [2 4]2((2 3)(gix -k 2±) + (3 4)A;2± v^^r) 

fc2±-S23 S34 5234^5?^ 

k3±k;^ {q2± k2±^/x^ - qi± {ql^ + kl^) ^/x^ + k2i_[2A]^/xZ) 

S2aa'S4bb'k2±\/x^ 

k2± {-<hi- <l2± /'L /'■2 ■I':', - - k2±)k:^^_ [3 l]/r7T7 + y2± ^2^ [2 1] .fi ) 

{k3±)'k;^ X2 ^ ^ |g2±P [24]x3 (-(23)^+ (34)7^) 

(34)/C2±S23 (|A:3±P X2 + |A;2^P Xg ) ^ (2 3)S34 S23i\/X4 

g2± [2 4]^ ((3 4) {ql^ + A;*^ ) ^ - (2 3)^2^ ^) 

S23 S34 5234^^2 a^4 

|gl±P (A:3±)^fe2± ^4± ^2 

^X3 X4 (34)A;2_lS23 (kl± - ?2±P + S234) (|^3±P X2 + |A;2±P 2^3 ) 
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X3^/x^{3 4:)k2±S23 S234 {\qi± - g2±P + S234) 
1 



, ^ ^ ,-=— (^3 {-X2 ih^ q2± [2 4] + [2 3] [3 4]) 

+ q2± fcsx k*^ [34]v/^r^-g2± |A;2±P [24]a;3 ) 

+ ql± h± X2 {q2± [2 4]x3 + /c3± [^MV^TW+h^ [2 4]x4 )) 



(6.14) 



/l33(gi,A;2+,A;3+,A;4,g2) = 2v^ 



''2± 



(gi± g2± (2 3)A;*^ 



_S2aa'(2 3)A;2±S34 y/x^ 

- (2 4)[3 4]7^(gi_L - k2±) (g2± + h± ) + k2^\/x^ {q2± S23 + (g2± + ^4± ) S34 )) 

ql^ iq2± + h± f {qi± K± - ^±['2 3]^^) |g2±P k3± y/X2 X3 



S2aa'S4bb'{'2 3)k2±y/x^ 



(2 3)fc2±S34 (1 - X2 



|gl±P (fc4±)' ((34)^3^ (24)fc,\ Ig^^p (gi^ -k2±)k;^X3 



a;4v/^(2 3)(3 4)A;2±S234 (|gi± - g2±P + S234) ' S2aa' ^2X^34 (1 - X2 

^ g4wfc4± v/^v/^+ (g2± + h± ) (((3 4)A:*^ + k^^ [3 4]) (2 4)fc;^ 

(2 3)^2X534 X4 

^ |g2±P [2 3] ((2 4)7^+ (3 4)7^) ^ |g2±P 7^^0:4 



(2 3) -534 S234 



(2 3)S34 {I-X2 



[2 3] 



(/^4±)^a;2 + g2± 



(2 3)A;2±S34 S234a;4 

- q*i± k2±h± ih± + q2± X4 )) 



(2 4) I A:: 



2± 



+ 



/X2 



{3 4)k2±k;^ 

\/xr 



(6.15) 



where in Eqs. ( |6.13|) -( |6TT5D we have used the three-particle invariants, S2aa' = (^2 — qiY 
and Sibb' = {h + q2Y- 

Eq. ( |6.11| ) must not diverge more rapidly than l/|q'j±| for 0, with i = 1,2, 

in order for the related cross section not to diverge more than logarithmically. Since 
Eq. (|6.11| ) is proportional to l/l^ixp, the NNLO Lipatov vertex must be at least linear 
in \qi±\, 

hm A'\q^,k^,\kl^^k^,^q2) = 0i\q.^\), (6.16) 
which is fulfilled by Eqs. (|6l2D-(|6J[5D. 



As a consistency check on Eq. ( 6.11 ), in the further limits 1/2 ^ 1/3 — or ?/2 — 1/3 ^ 
?/4, the NNLO Lipatov vertex in Eq. ( |6.11| ) must factorize into a NLO Lipatov vertex 
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convoluted with a multi-Regge ladder, 



0-G53 



(6.17) 



with gi2 = Qi — k2, and 



1 



(6.18) 



1 



,gfc"d,c- C3^qu,k^,\q2) 



with gi2 = ^2 + ^4- 

In the triple collinear limit, /c2 = Z2P, k^ = z^P and k^ = z^P, with 22 + -2^3 + -24 = 1, 
the coefficients of the NNLO Lipatov vertex ( |6.12| )-( prT5| ) reduce to the splitting functions 
( |5.52| ), and amplitude (|6.11|) factorizes into a multi-Regge amplitude ( |4.11| ) times a double- 
coUinear factor ( p.49[ ) 

lim A^^-^^{p:-X' I k^2\K\ kT I p7,Pt) = 
K2IIK3IIK4 

yA^'-''{p:% k\^ I p'' I piry,^) ■ spiit^-:''{k^,\k'2\ ^3^) • 



6.4 NNNLO impact factors in the quadruple collinear limit 

In the quadruple collinear limit, ki = ZiP, with zi + Z2 + ^3 + 2:4 = 1 a generic amplitude 
is expected to factorize as 



lim j\_-did2d3d4 

ki\\k2\\k3\\k4 



M 5 ■■■) 



11J2 I.U3 11 
■■■,1^1 ,1^2 ' ""3 ' "'4 

yA-'-{..., P\ ...) ■ Split{7^^^^^^^^(^r, A;2^% A;3^% ^I^) • 



(6.19) 



Accordingly, we show that we can write Eq. (|6.8|) as 

lim ^^^-^^ (p^'' , A;^ , ^2^^ , k'^^ , A;^ | , ) 

A:i||fc2||fc3||fc4 



(6.20) 



by taking the quadruple collinear limit of the NNNLO impact factor. 

In the quadruple collinear limit, the functions A^'^^ of Eq. (|6.5| ) yield a cubic divergence 



as S1234 = {ki + /c2 + ^3 + ^4)^ ^ or Sijk 0, or Sij —>■ with i,j,k = 1,2,3,4. 
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Analogously to Sect. |5.(j| , a function A^'^s differs from its reflection by a term which 
contains only a quadratic divergence in the vanishing invariants. Using this property 
and Eqs. ( p. 61 ) and (|6.7D, we obtain a reflection identity and dual Ward identities, up to 
quadratically divergent terms, 

/1S;49('„i^a. 1,1^1 1.1^2 Lf3 l.'^4\ __ Ag-Agf^fa. 1.1^4 1.^3 l.'^2 (P,91\ 

A9'A9(r,'^a. 1.^2 U^i t.l^4\ _|_ AdAgf^i'a. Ufl t.1^2 t.1^4 UKi\ 

^ \Pa 5 '^1 ? '^2 ' '^3 ' '^4 J ^ iFa ) '^1 5 '^2 ' '^4 ' '^3 / 

+A^-''%p:-, k\\kl\ kl\ kf) + A^-^'^{p':-, kl\ k^,\k^,\ kl^) = 0, (6.22) 

^ \Pa 1 '^1 ) '^2 ' '^3 ' '^4 / ^ \Pa i '^1 ; '^3 ' '^2 ? '^4 / 

\Pa 1 '^Z 1 '^1 1 '^2 1 '^A ) ^ ^ \Pa 1 '^1 5 '^3 ' '^4 5 '^2 J 

fc3^^ fcr, C A:^^^) + A^-^''{p7: kl\ kl\ k'i\k'^-) = . (6.23) 

We note, however that the last identity is not independent from the first two. Using the 
identities ( |6.21|) -( |6r23D in Eq. (|6.8| ), we can factorize the color structure on a leg 



(Tes'4 

= ^9r'" 1/$^ (F'^-F'^-F'^-3),,^A^^'^(p:^fc:r,C%C%0 |>(6.24) 
thus amplitude (|6.^) can be put in the form of Eq. (|6.2U|) with coUinear factor 



Split?.j^ ^ {ki^ , k2^ , k^^ , k^*) 

= g'J2 (F'^^F'^^^F'-s)^^^ sp\it^-:'%kZ\kZ\kZ\kT) . (6.25) 



The splitting factors splitl^^^ are the functions A of Section |6.1| in the quadruple collinear 
limit, up to quadratically divergent terms, and thus they fulfill the identities, Eqs. (|6.21 



(|6.23| ). The splitting factors of PT type can be soon read off from Eqs. (|6.3|) - (|6.5|) , while 



the ones of non-PT type can be given in terms of three functions of the collinear momenta. 



splitr"^'^(A;+,fc^-,fc+,A;+) = 2^2—4 



1 



yilil(12)(2 3)(3 4) 

^2 I 

spUr''{k^,kt,kt,kt) = 2v^^ (12) (2 3) (3 4) 
splitr''(A:i+,/^2",/^3'",A;+) = 2v^y^ (1 2) (2 3)(3 4) 
spUi-^'' {kt,kt,k^,kt) = 2v^^(i2)(23)(34) 

split'T^' (k+kt,kt,k7) = 2^2^^——^—— (6.26) 
^ + ^ ^' ^' ^' yirzl(12)(2 3)(3 4) ^ ^ 
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o^il L _ V I 


A-+ 

"-2 ) 




K4 ) — 
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o^il t V 1 


1^2 1 
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-^2(4,3,2,1) 




Spllt_ (^/C]^ 


^+ 

1^2 1 




K4 ) - 


^3(1,2,3,4) 




spiitr^^(A;r 






K) = 


-^3(1, 2, 4, 3) + B2(3, 2, 4, 1) + -62(3, 2, 1, 4) 


splitr^^(A;+ 


^2 ) 


^3 


kt) = 


-B2(l,4,2,3) + B3(3,4,2,l) 


--62(4,1,2,3) 


splitr^^(A;+ 


^2 ) 


^3" 




-^3(4,3,2,1) 




splitr^^(A;+ 


^2") 


^3 


K) = 


^2(1, 2, 3, 4) 





with 



-61(1,2,3,4). 
(12) (2 3) S34 



-235(1,2,3) 



V^2^3^4(12) 



+ 



^3 



V^(^3 + ^4)[12] (1 - Zi) (^3 + 2^4) V ^1^4 

e(l,2,3,4)'((24)[12] + (34)[13])[34] , /i^6(l, 2, 3, 4)'[2 3] 



(1-^4)6(1,2,3,4) 



+ 

•51234S234 V ^1 

v/il(l 2) [2 3] ( Vir(l 4) - (2 - z,) 6(1, 2, 3, 4)) 

(1 - 2;i)s234 



S234 



(6.27) 



Z2^/^ ^^izi + Z2)S{l,2,3) 5(l,2,3)'e(l,2,3,4)[34] 



-62(1,2,3,4) 

2^/2 

[34]S23 (l-^l)V^ V^l^Sl2 Si234V^'5l2 S123 

(5(1, 2, 3) (v/il(5(3, 4, 1) + ^5(3, 4, 2)) + ./z^z^s^ su 

■Sl234v^(12)[3 4] 

, ;2iV^5*(3,4,2) ^i(34)[12]((13)[23] + (14)[24]) 



(1 - 2;i)s234 



S1234S234 



(6.28) 



-63(1,2,3,4) 



2V2 



Z2\fZ2Zi 



Z2 



S\2 S23 S34 

(Vi^(12)-^4(13))5(2,4,3) 
2 



2:1(23)534 ^(12)S23 



^1(2:1+ 2:2) V 1-2:4 2:3 + 2:4 

2:2(75^6(1,2,4,3)523 - (2 3) (3 4) [2 4]) 



Z\Z4 (2:1+2:2) 

(^(1 2)5(2, 4, 3)' + v/ir(l 3) (^2 (1 2) (2 3) + ^^(1 3) (3 4) + ^4(1 4) (3 4)) [2 4]') 



+ 4/^2:2(13)5234 



^1*1234 



Z2S23 S34 5*{1, 3, 2) , 512 5(2, 4, 3)^6(2, 3, 4, 1) 



Zl (1 - 2:4) [2 3]5i23 



+ 



•5l234V^l'^234 



•5l234'5l23 



2:4(1 3)'[2 4] (6(1, 2, 4, 3) [2 3] - 6(2, 3, 4, 1)[1 2]) 
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^(13)^ (-v/i^e*(2, 3, 4, 1) + yilil[24]) 

Zl Si23 

with (5(1, 2, 3) as in Eq. (|5.51|) , and 



(6.29) 



6(1, 2, 3, 4) = v^(l 4) + v/i^(2 4) + ^^{'i 4) . (6.30) 

As in Section p.6| , summing over the hehcities of gluons 1, 2, 3 and 4, one can obtain 
the two-dimensional polarization matrix, 

spiit^j:^^(fcr, kf, A;r)[spiit?_-:''(A;r, ^r, ^d]* = 5^^'4^PaV , (6.31) 

where P^'^^ = P^^^ , and P^^^ = (P^I*^^)*. Averaging then over the trace of matrix 
( |6.31| ), i.e. over color and helicity of the parent gluon, one can obtain the unpolarized 



Altarelli-Parisi gluon triple-splitting function 

^ ^ \^^\\ir,'\k'i\kl\kf,kf)\^ = ^{P^-'<^), (6.32) 

^ ' UU1U2U3U4 -L^J-i 



2(iV2 



with (ps^^s) = P^'^^ . As in Section |5.6| , the sum over colors can be done using Eq. ( p.5| 
and we obtain. 



>'(T4 ; I , 



|Split3-^^^^^394(^..^ ^..^ ^.3^ ^.4)|2 ^ ^ |split^_7'^(CS C% C% 

o"£54 

(6.33) 

with C^{Nc) as in Eq. ( p.8| ). It is then clear that for the splitting functions ps^"^, with 
n > 4, the color will not factorize since LCA, Eq. (|2.5|) , is not exact any more. We do not 



compute here P^^^ and P_^r^^, all the information about them being already contained 
in Eqs. dOeQ-dOOl). 



7 Conclusions 

In this paper, the structure of QCD amplitudes in the high-energy limit and in the coUinear 
limit has been explored beyond NLO. We have computed forward clusters of three partons 
and four gluons, which in the BFKL theory constitute the tree parts of NNLO and NNNLO 
impact factors for jet production. In the BFKL theory the NNLO impact factors could 



be used to compute jet rates at NNLL accuracy. In Sect. pTll , pl2| , |5l3| and pl^ , we have 
computed the tree parts of the NNLO impact factors for all the parton flavors. On these 
we have performed in Sect. a set of consistency checks in the high-energy limit, and we 
have obtained in the triple collinear limit (Sect. ^.6|) the polarized, the spin-correlated and 
the unpolarized double-splitting functions. The last two agree with previous calculations 
by Catani-Grazzini and Campbell-Glover, respectively. They can be used to set up general 
algorithms to compute jet rates at NNLO. 
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From the four-gluon forward cluster we have obtained in Sect. |6.1j the tree part of 
the purely gluonic NNNLO impact factor. In the quadruple coUinear limit, this yields 
(Sect. the purely gluonic unpolarized triple-splitting functions. They could be used 
to compute the three-loop Altarelli-Parisi evolution, or to compute jet rates at NNNLO. 
In addition, by separating a central cluster of three gluons out of the four-gluon forward 
cluster, we have computed the emission of three gluons along the ladder, Eqs. ( |6.11|) - (|6.15|) , 
which contributes to the NNLO Lipatov vertex. This constitutes one of the universal 
building blocks in an eventual construction of a BFKL resummation at NNLL accuracy. 

Finally, inspired by the color structure in the high-energy limit, we have found a 
compact color decomposition of the tree multigluon amplitudes in terms of the linearly 
independent subamplitudes only, Eq. (|2.9D. It would be interesting to analyse whether 
this structure generalizes to multigluon amplitudes at one loop, and beyond. 

The decomposition in rapidity of amplitudes in terms of gauge-invariant parton clus- 
ters performed in this work suggests naturally a modular decomposition of a generic 
multiparton amplitude, where each module is an n-parton cluster. Such an approxima- 
tion could be tested against existing approximations of multiparton amplitudes ^ 



In the high-energy limit, the cluster decomposition seems superior, in that it does not use 
only PT-type subamplitudes, like the Kunszt-Stirling approximation [^, and within a 
cluster it is not limited to collinear kinematics, like the Maxwell approximation [R^]. 



A Multiparton kinematics 

We consider the production of n partons of momentum with z = 1, n and n > 2, in 
the scattering between two partons of momenta Pa and p^. 

Using light-cone coordinates p"^ = po ± Pz, and complex transverse coordinates p± = 
px _|_ j^pp^ with scalar product 2p ■ q = p^q^ + — p±q*l — p*j_q±, the 4-momenta are. 

Pa = (p+/2,0,0,p+/2) = (p+,0;0,0) , 

Pb = (p,-/2,0,0,-p,-/2) = (0,p,-;0,0) , (A.l) 

Pi = {{pf +Pi)/'^^ Re [pi±] , Im [pij_] , {pf -p~)/ 2) 

= {\Pi±\e^',\Pi±\e~'^";\Pi±\cos(j)i,\pi±\sm(j)i) , 



where y is the rapidity. The first notation in Eq. ( |A.1[ ) is the standard representation 



p/j. _ {jp ^p^ ^py ^p^)^ while in the second we have the -|- and - components on the left 

^By convention we consider the scattering in the unphysical region where ah momenta are taken as 
outgoing, and then we anahtically continue to the physical region where < and p° < 0. Thus partons 
are ingoing or outgoing depending on the sign of their energy. Since the hehcity of a positive-energy 
(negative-energy) massless spinor has the same (opposite) sign as its chiraUty, the heUcities assigned 
to the partons depend on whether they are incoming or outgoing. Our convention is to label outgoing 
(positive-energy) particles with their hclicity; so if they are incoming the actual helicity and charge is 
reversed. 
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of the semicolon, and on the right the transverse components. In the following, if not 
differently stated, pi and pj are always understood for 1 < i, j < n. From the momentum 
conservation, 

n 

= ^Pi±^ 

i=l 

n 

Pa = -Y.pt ^ 

i=l 

n 

Pb = - J2pi ^ 

i=l 

the Mandelstam invariants may be written as, 

Sij = 2pi ■ Pj = pfpj + prp+ - pi^p*^ - p*^_Pj± . 

so that 



S = '2Pa-Pb=Yl PiPj 

n 

•ai = 2p„ = (A.3) 



n 

Sbi = 2pb-pi^ - Y PiPj ■ 



Massless Dirac spinors ^±(p) of fixed helicity are defined by the projection. 



Mp)-^^Hp), (A.4) 



with the shorthand notation. 



Mp) = b±), ^±(p) = (p± I , 

(pk) = {p - \k+) = ^_{p)4,+ {k) , (A.5) 
\pk] = ip+ \k-) = i)+{p)i)_{k) . 

Using the chiral representation of the 7-matrices, 
and the normalization condition: 

(p±|7^|p±) = 2p^, (A.7) 
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and the complex notation = Ipxle**^, the spinors for the momenta ( |A.1| ) are|^ 



'Pi e 




i(t>i 



1p+{Pa) = i 



\ / 





v 








ip-{Pa) = i 



J 



( \ 



\ -VW ) 

( \ 







(A.8) 



/ \ 



^~{pb) = -I 



( \ 





Ph 

\ / 



V / 

Using the above spinor representation, the spinor products for the momenta ( |A. 1| ) are 



{PiPj) = Pi± 





'til^Pi. 
Pi 



(A.9) 



-Pb Pi 



{PaPi) 

{PiPb) 
{PaPb) 

where we have used the mass-shell condition \pi±\'^ = P^P^ ■ Note that in the present 
convention the spinors ( [A .81 ) and the spinor products ( |A.9| ) differ by phases with respect 
to the same in Ref. 



We consider also the spinor products {pi + ■ Pk\Pj+) i which in the spinor representation 
( |A.(j| ) take the form, 



{Pi + \l-Pk\Pj+) 

{Pi + h-Pj\Pa+) 
{Pi + \l-Pj\Pb+) 



{ptptPk - PtPj±Pl± - Pi±Pj'Pk± + PiLPj±Pt) . V/C 



pTpJ 

-Pi 
pf 



{pfPj -Pi±Pj±) , Vj 



(A.IO) 



-Pb 



Pi 



-pTPj± + Pi±Pj) 



"The spinors of the incoming partons must be continued to negative energy after the complex conju- 
gation. For instance, = i (^\/~-pf, 0, 0, 0^ . 
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The spinor products fulfill the identities {i = Pi,j = pj] 





= -m 




= -m 




= sign(p°p°) [ji] 


{{^ + \Y\J+)r 






2pj ■ Pj Sij 




= m{kj) 




= mm 




= {tk){jl) + {U){kj) 


m m 


= mm + mikj] 



(A.ll) 



and if X^iLi = then 

n 

Y,mm = o. (A.12) 



1=1 



Throughout the paper the following representation for the gluon polarization is used, 
which enjoys the properties 

et*{p,k) = e;{p,k), 

et{p,k)-p = e^{p,k)-k = 0, (A.14) 

u=± P 

where k is an arbitrary light-like momentum. The sum in Eq. ( |A.14| ) is equivalent to use 
an axial, or physical, gauge. 



B Multi-Regge kinematics 

In the multi-Regge kinematics, we require that the gluons are strongly ordered in rapidity 
and have comparable transverse momentum, 

>...>?/„; ~ ... ^ b„_L| . (B.l) 

Momentum conservation ( [A.2| ) then becomes 

n 

= X^m, 

i=l 

Pt ^ -Pt, (B.2) 

Pb - -Pn ■ 
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The Mandelstam invariants ( [A .31 ) are reduced to 



S = 2pa-pb^ ptPn 
Sai = '^Pa ■ Pi ^ -ptPi (B.3) 
Sbi = 2pb ■ Pi ~ -ptPn 



Sij = 2pi-pj ~ \Pi±\\Pj±\e 
to leading accuracy. The spinor products ( |A.9| ) become 



IVi-Vjl 




{piPj ) ~ -\l-^Pj± ioTyi> yj 

(PaP.) ~ (B.4) 



p. 



{PiPb) ^ i\/ptPn . 
(PaPfe) ^ -\/ptPn ■ 

C NLO Multi-Regge kinematics 

We consider the production of n partons of momenta pi, ...jPn, with partons 1 and 2 in 
the forward-rapidity region of parton pa, 

~ ?/2 > 1/3 > ••• > Vn ; Iml - Iml - ••• - \Pn±\ ■ (C.i) 

Momentum conservation ([A.2|) becomes 

n 

= X^m, 

i=l 

p+ ^ -ipt+pt), (C.2) 

Pfc" - -Pn- 

The spinor products ( |A.9| ) become 



(PaPfe) = -V^ ^ -V (P^ + Pt)P 



(PaPn) = — ^ Pn± ^ iP^APaPb) 
Pt \Pn±\ 



Pk V Pfe 



(PfePb) = H -Pb Pk - h PkPn > /C = 1, . . . , n - 1 
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{VnVb) = isJ-VbVt -AVnl. 
{PkPn) = Pkl. 




-Pr. 



{P1P2) = Pl±\ —r-p2± 




(PkPi) = Pk± 






Pi 




(C.3) 



-Pn 




k = 1, . . . ,n — 1 



-Pii 




Pi 
+ ' 



/c=l,2;z = 3,...,n — 1 



which differ by phases with respect to the same spinor products in Ref. because of 
the convention for the spinor representation we use in Sect. 0. 



D NNLO Multi-Regge kinematics 

The extension to the production of n partons of momenta pi, with partons 1, 2 and 

3 in the forward-rapidity region of parton pa-, 

?/l ~ ?/2 ^ 1/3 > Z/4 > ••• > Vn ; \Pll\ ^ Iml - ••• - \Pn±\ , (D.l) 

is straightforward. We mention it here because by taking the further hmit ?/i ^ ?/2 — Z/s, 
one obtains the kinematics of the NLO Lipatov vertex (sect. |5.5| ). 

With Eq. ( p.l| ), momentum conservation ( [A.2| ) becomes 

i=l 

pI - -{Pt+Pt+Pt)^ (D.2) 

Ph - -Pn- 

The spinor products (|A.9|) become 



{PaPb) = -\fs ^ -\ {pt +pt + Pt)Pn 



{PaPn) = -i\I^^Pn±^i-r^{PaPb) 
PZ \Pn±\ 



{PaPk) = -i\I^T~Pk±--Hl ^ ^+ Pfc± , A; = l,...,n-1 

Pk V Pk 



{pkPb) = iy-PbPk -'^yPkPn^ k = l,...,n-l 

{PnPb) = i\[-PbPt -APnA^ 

{PkPn) = Pfc± -Pn±W^ ^ -Pn± , /c=l,...,n-l 

\ Pk \ Pn \ Pt 
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{VkVi) = 




/c = l,2,3;i = 4, ...,n — 1, 



while the others spinor products remain unchanged. The spinor products ( p. 3D generahze 
straightforwardly to the kinematics ( |6.1|) . 



E The Sudakov parametrization 

We want to elucidate the relationship between our parametrization of the momenta and 
the one of Ref. [0. Recalling the last of Eqs. ([A.l|) , we can write, 



Pi 



^(1, 0, 0, 1) + (0, Re[p,J, Im[p,J, 0) + ^^(1, 0, 0, -1) 



(E.l) 



where is the sum of the three momenta, the Xi are the momentum fractions and we 
used the mass-shell condition plp~ = \pi±\'^- This is exactly what is obtained from the 
general Sudakov parametrization of Ref. ||11||, 



Pi 



XiP" + k^\^ - 



Xi 2p ■ n 

through the following choices for the lightlike vectors. 



and the identification. 



;i, 0,0,1) 



and 



1,0,0,-1) 



fc^^, = (0,Re[piJ,Im[m],0). 



(E.2) 



(E.3) 



(E.4) 



The spin-correlated splitting functions of Ref. []TT1 are expressed in terms of the vectors ki 
defined as k^ = k'^^ — ZiP'^, where, as in our case, the Zi variables represent the momentum 
fractions in the coUinear limit. In order to compare Eq. (|5.59|) with the spin-correlated 



splitting functions of Ref. |11], we must project the latter onto the helicity basis, namely 
to contract them with the polarization vector, Eq. ( |A.13| ), 



e±(P,n) = -^(0,1,^^0) 



(E.5) 



The contraction of the fcf vectors with is. 



ki ■ 



(E.6) 



with i, j, / = 1, 2, 3 and j, / ^ i, with the analogous expressions for e obtained by complex 
conjugation. 

For the off-diagonal terms, P3^9f2h^ g^^j relative minus sign between the results 



of Ref. |]TT| and ours, which, however, has no physical relevance. 
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